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1.  INTRODUCTION 

The  goal  of  the  present  work  is  the  integration  of  active  and  passive  control  of  structural 
vibrations.  First  a  procedure  is  developed  to  achieve  maximum  passive  viscous  damping  in  the 
frequency  range  of  interest.  Then  the  active  control  is  added  to  increase  the  effective  damping  to 
the  level  required. 

Passive  damping  is  either  viscous  damping  or  that  obtained  with  viscoelastic  materials.  It 
is  representable  by  devices  connected  between  two  structural  nodes  or  between  one  node  and 
ground.  It  is  achieved  mechanically  by  special  truss  members  which  incorporate  dashpots, 
electo-mechanical  devices,  or  viscoelastic  elements  where  the  viscoelastic  material  is  deformed 
in  shear.  Viscous  damping  is  represented  in  matrix  form  as  F  =  C  z  .  Viscoelastic  damping  is 
represented  herein  by  (1)  using  a  complex  Young’s  modulus  which  is  frequency  dependent 
[N3],  where  the  F  =  G(0))  or  (2)  by  using  the  mini-oscillator  approach  to  represent  the 
frequency  dependence  of  the  complex  Young’s  modulus,  as  presented  by  McTavish  and  Hughes 
[M2  ].  The  fractional  calculus  method  propounded  by  Bagley  and  Torvik  [B1  ]  is  of  interest,  but 
has  not  been  pursued  in  the  present  study.  The  active  control  investigated  herein  is  mainly 
constant  gain,  closed  loop  control.  It  is  assumed  to  be  achievable  by  control  forces,  in  matrix 
form,  of  F  =-  G  z  •  H  z  . 

The  second  chapter  presents  analytical  background.  Much  of  control  theory  is  written  in 
state  space  form,  which  for  structural  problems  means  reformulating  n  second-order  differential 
equation^  into  2n  equivalent  first-order  equations.  Part  of  the  cost  is  that  n  redundant  first-order 
equations  have  to  be  introduced  and  a  corresponding  amount  of  redundant  information  can  be 
produced.  So  some  of  the  features  of  eigenproblems  for  both  second-order  and  state  space 
differential  equations  are  reviewed.  Also,  important  relationships  are  presented  between  the 
upper  and  lower  portions  of  the  right  and  left  eigenvectors  which,  though  easy  to  derive,  are  not 
readily  available  in  the  literature. 

A  great  deal  of  emphasis  is  put  on  sensitivity  of  the  eigenproperties  to  damping  and 
control  forces.  This  is  summarized  in  the  third  chapter.  The  equations,  initially  derived  by 
others,  for  calculation  of  eigenvalue  and  eigenvector  sensitivities  are  summarized  including  some 
examples.  A  unique  feature  is  the  introduction  of  sensitivities  calculated  also  from  the  invariants 
associated  with  the  similarity  transformation  between  system  coordinates  and  modal  coordinates. 

In  the  fourth  chapter,  optimization  of  passive  damping  is  summarized.  A  ten  bar  truss  is 
used  as  a  practical  example  and  the  optimization  of  both  viscous  damping  and  viscoelastic 
damping  is  discussed  in  detail. 
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Chapter  five  deals  with  active  damping  and  control  and  design  of  the  [C+G]  matrix  to 
achieve  the  desired  combination  of  passive  and  active  control.  The  existing  theory  of  natural 
modal  control  is  reviewed.  Two  approaches  are  presented  for  designing  the  combination  of 
passive  viscous  damping  and  full  state  feedback  control.  The  first  is  iterative  and  makes  use  of 
eigenvalue  sensitivities.  The  second  is  based  on  a  perturbation  formulation,  where  the  desired 
changes  in  eigenvalues  are  specified  and  the  corresponding  changes  in  the  closed  loop 
characteristic  matrix  in  state  space  are  predicted.  Restrictions  on  the  choice  of  eigenvalues  and 
eigenvectors  are  discussed  in  light  of  related  literature  on  eigenstructure  assignment.  The  concept 
of  cost  functions  is  reviewed  and  the  relationship  of  solutions  of  the  Riccati  equation  to  the 
present  work  is  outlined. 

2.  ANALYTICAL  BACKGROUND  FOR  MODAL  THEORY 

The  purpose  of  this  chapter  is  to  review  the  modal  theory  and  notation  for  structural 
vibration  problems,  with  emphasis  on  assumptions  regarding  symmetry  and  proportionality  of 
damping  or  control.  First  the  second  order  differential  equations  of  motion  are  outlined;  second, 
the  state  space  form  of  the  equations. 

2.1  Second  order  differential  equations  of  motion. 

2.1.1  Eigen  properties  and  orthogonality  with  no  damping. 

With  regard  to  nomenclature,  a  bold  capital  letter  or  any  symbols  in  brackets  [  ] 
represent  a  rectangular  matrix.  A  vector,  or  column  matrix,  is  represented  by  a  bold,  lower  case 
letter  such  as  z,  symbols  in  braces  {  },  or  an  underlined  Greek  letter,  as  with  A  diagonal, 

square,  matrix  is  represented  by  a  superscript  "D"  on  brackets,  that  is  [ 

The  differential  equations  of  motion  for  the  case  of  no  damping  are 

Mz  +  Ki=F(t)  (2.1) 

Mz  +  K  z  =  -Gz  -Hz.  (2.2) 

where  M  is  positive  definite  and  symmetric,  K  is  indefinite  and  synunetric  and  F(t)  =-  Gz  -Hz. 

1  lie  G  and  H  are  control  matrices.  The  eigenvalue  problem  is 

[-Ma)J+Kl(e^  =  0  (2.3) 

The  eigenvectors  may  be  assembled  into  an  n  x  n  <I>  matrix.  The  orthogonality  relationships 
are 

<d’^M<I>  =  (2.4) 
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and 


(2.5) 


<I)^K<1>=(M  0)^1 
‘  n  n' 

2.1.2  Reconstruction  of  the  M  and  K  matrices  from  eigenvalues  and  eigenvectors. 


Note  that  the  M  and  K  matrices  may  be  reconstructed  by  using 


(2.6) 

and 

(2.7) 

Further,  we  see  that 

(2.8) 

and 

m  ‘k  =  <I>  [0)^1  ^d>'' 

(2.9) 

2.1  J  Uncoupling  of  the  differential  equations  and  modal  control. 

The  differential  equations,  Eqs.  (2.1)  and  (2.2),  are  uncoupled  by  the  coordinate 
transformation 

z(t)  =  <I>q(t)  (2.10) 

T 

and  pre-multiplying  Eq.  (2.1)  by  d>  . 

<I>‘^M<h  q  +  q  F(t)  =  N(t)  (2. 1 1 ) 

IM  q  +  [Mf  =  -  [g]  q  -  [h]  q  =  N(t)  (2.12) 

The  diagonal  matrices  on  the  left  side  of  Eqs.  (2.12)  indicate  decoupling,  but  if  [g]  and 
[h]  are  non-diagonal,  the  equations  are  still  coupled.  If  [g]  and  [h]  are  diagonal,  which  is  called 
natural  control  [M3],  then  the  equations  are  uncoupled. 

If  the  control  matrices  in  N(t)  are  specified,  and  O  F(t)  =  N(t),  then  F(;)  is  found.  It  is 
efficient  to  use  the  orthogonality  relationship,  Eq.  (2.4),  to  avoid  inverting  <1>  . 

F(t)  =  N{t)  =  MO)  (M  N(t)  =  -  MO  [M  [[  g]  q  +[h]  q]  (2.13) 
?hanging  back  to  the  z(t)  coordinates 

F(t)  =  N(t)  =  MOl  N(t)  =  -M  <DI  [[  g]  d>**  z  +[hl  O'^z]  (2.14) 
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The  differential  equation  with  feedback  is  now 


Mz  +  M  d>[  M [[  g]  d>  '  i  +(h]  d>  *z]  +  Kz  =  0. 

(2.15) 

which  is 

Mz  +  Gz  +[H  +  K]  z  =  0. 

(2.16) 

where  G 

=  M  d)[  MJ‘^[  g]  d>‘‘  and  H=  M  d)[  h]  d>  * 

(2.17) 

Pre-multiplying  Eq.  (2.16)  by  M** 

z  +  d)(  [[  g]  d>'‘  z  +[h]  d>  *z]  +M''  Kz  =  0. 

(2.18) 

or 

z  +  M‘‘  Gz  +M''lH  +  KJ  z  =0. 

(2.19) 

which  is  of  the  form 

I  z  +  Cj  z  +  KjZ  =0. 

(2.20) 

2.1.4  Criterion  for  proportional  damping  or  proportional  control. 

The  question  now  addressed  is  whether  the  matrices  H  and  G  in  Eq.  (2.19)  represent 
proportional  control.  The  test  which  applies  is 

9 

is  Cj  Kj  =  KjCj  (2.21a) 

m  ‘gM'‘  [H  +  K1  =  M‘^  (H+Kl  M'‘ G  (2.21b) 

Since  G  and  H  can  be  chosen  separately,  the  Eq.  (2.21b)  is  separated  into  two  parts: 

9 

M‘‘g  M'‘  H  =  M*^  H  M'‘  G  (2.22a) 

9 

M'‘g  M'‘k  =  M'*K  M'*  G  (2.22b) 

Now  substitute  for  !VI'*G  and  M  *  K  from  Eqs.  (2.18)  and  (2.19) 

9 

<1>[  g]  4)'‘  d)  I  ci);|  =  cp  [  (D^i  M'*  M  d>[  MJ'^[  g]  d)'*  (2.23a) 

which  reduces  to 

[  [  g]  [  0)^1  ^  =  1 0)2]  ^  [  g]  (2.23b) 

The  answer  is  that  Eq.  (2.23b)  does  not  represent  an  equality.  However,  the  two  sides  of 
the  equation  will  be  equal  if  [gj  is  diagonal,  which  is  then  a  requirement  for  proportional  control. 
Then  Eq.  (2.23b)  becomes  the  following  equality 

I  Mj  ’’  [  gl'’  [  <0^1  “  =  I  «>^l  “  I  M„l'°  I  8l°  (2-24) 
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The  advantage  of  proportional  damping  or  proportional  control  is  that  the  eigenvectors  <I>  do  not 
change,  which  means  that  the  coordinate  transformations  which  decouple  the  original  matrices 

will  also  decouple  the  closed  loop  equations.  A  direct  way  to  check  that  this  is  true  is  to  simply 

T  T 

look  at  the  products  <1>  G  <!>  and  d>  H  <I>. 

The  latter  product  <I>'^H<1)=  M  <1>[  h]  <I>’  *  cD  (2.25) 

reduces  to  (hj ,  (2.26a) 

and,  similarly,  <I>^G<I>=  [gj.  (2.26b) 

So,  the  open  loop  eigenvectors  decouple  the  closed  loop  equations,  but  only  if  [h]  and  [g]  are 
diagonal  matrices. 


2.2  State  space  form  of  differential  equations. 

2.2.1  Eigen  properties  for  symmetric  plant  and  control  matrices. 

The  differential  equations  in  state  space  form  are 

M*  q +K*  5  =  -  G*  n 
I  n +A^q  =  .[M*r' 

which  defines  A  .  With  G*  =  0,  Eq.  (2.27)  is  written  in  detail  as 
0 

'  M  0  1  [  z  (t)  1  r  C  k1  (  z  (t)l  0 

-  >4-  -  > 

_0  -KJ  lz(t)J  [K  oj  lz(t)J  lo. 

The  open-loop  eigenproblem  is 


(2.27) 

(2.28) 

(2.29) 


[IX-  A^  ]  \|/  =0  (2.30) 

Xt 

and  the  solution  to  the  homogeneous  equations  is  q  ^  e  *■  ,  where  is  the  xth  eigenvector. 

The  2n  eigenvectors  may  be  arranged  in  columns  in  a  2n  x  2n  matrix,  designated  T  ,  and  the 
following  orthogonality  relationships  [N3]  exist 


y  =  'P^K*H'  =  -B^A  (2.31) 

where  the  eigenvectors  may  be  normalized  so  B*^  =I .  The  matrix  A  is  a  diagonal  matrix  of 
eigenvalues,  X^. 


The  y^  are  sometimes  called  the  right  eigenvectors,  because  many  workers  in  the  field  of 

control  seem  to  be  in  the  habit  of  using  right  and  left  eigenvectors,  even  though  left  eigenvectors 
are  not  needed  for  decoupling  if  the  plant  and  control  matrices  are  symmetric.  The  concept  of 
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left  and  right  eigenvectors  is  explained  in  the  next  section. 

2.2.2  Right  and  left  eigenvectors 

When  both  right  and  left  eigenvectors  are  used,  the  convention  herein  is  to  designate  the 
xth  right  eigenvector  by  and  the  ith  left  eigenvector  by  .  The  left  eigenvectors  arc  the 
eigenvectors  associated  with  the  transpose  of  the  A  matrix. 

[IX-A^]v=0  (2.32) 

The  right  and  left  vectors  are  arranged  in  columns  to  form  the  U  and  V  matrices,  of  size  2n  x  2n. 
Their  properties  are  such  that 

A^  U  =  A  (2.33a) 

and  U  is  diagonal.  (2.33b) 

Two  normalizations  are  usually  performed.  The  first  involves  the  right  vectors.  There  are 
various  choices  for  normalizing.  If  M  and  K  are  symmetric,  11  =  4^  and  the  right  eigenvec¬ 
tors  may  be  normalized  using,  from  Eq.  (2.31), 

u'^M*U  =  b‘^  =  I.  (2.34) 

In  a  more  general  situation,  if  U  M  U  is  not  diagonal,  the  individual  vectors  may  be 
normalized  by 

u  M  u  =1  (2.35) 

r  a  r  '  ' 

* 

where  M  is  an  arbitrary  matrix,  to  be  defined.  Among  the  possibilities  are  M  =  M  or  M  =1. 

The  lower  half  of  the  U  matrix  represents  the  displacement  mode  shapes.  In  Chapters  3 
and  5  a  special  normalization  involving  only  the  lower  half  of  the  U  matrix  is  found  to  be 
advantageous. 

Having  normalized  the  right  eigenvectors,  the  left  are  usually  normalized  by  making  the 
diagonal  matrix  indicated  in  Eq.  (2.33b)  equal  to  the  unit  matrix,  or 

U  =  I  (2.36) 
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2.23  Relationships  between  upper  and  lower  eigenvectors. 

The  complete  eigenvector  matrices,  U  and  V,  in  the  state  space  formulation  of  structural 
dynamics  problems  are  of  size  2n  x  2n  while  those  for  the  second  order  differential  equation 
form  are  only  n  x  n.  Further,  the  eigenvectors  in  state  space  occur  in  complex  conjugate  pairs 
for  underdamped  systems.  The  relationships  between  the  various  forms  of  the  vectors  has  not 
been  documented  well  in  the  literature.  The  purpose  of  this  section  is  to  outline  some  of  the 
relationships  between  the  portions  of  the  eigenvector  matrices.  The  right  eigenvector  matrix ,  U, 
is  first  divided  into  two  n  x  2n  matrices  designated  and  ,  for  the  upper  and  lower  portions 


of  the  U  matrix. 


U  = 


u 

u 

L 


(2.37) 


From  the  characteristic  equation  with  A  denoting  a  general  characteristic  matrix. 


A  U  =  U  A  (2.38) 

and  filling  in  the  submatrices  of  the  A  and  U  matrices,  the  following  two  sets  of  equations  result 

-C,  U  -K.U  =U  A  (2.39a) 

and  A  .  (2.39b) 

The  second  equation,  (2.39b),  shows  how  simply  the  upper  right  vectors  are  related  to  the  lower 
right  vectors  and  is  true  even  if  Cj  and  are  unsymmetric. 

The  left  eigenvector  matrix,  V,  is  subdivided  similarly  and  substituted  into  the 
characteristic  equation 

a'‘'v  =  VA  or  v’*'a  =  AV’^  (2.40) 

to  yield  =  ^'^u 

and  -Vj  K,  =  A  V^.  (2.41b) 

Here  the  relatidnship  between  the  upper  and  lower  left  vectors  is  not  so  simple,  unless  Cj  =0, 
which  is  a  situation  discussed  below,  under  section  2.2.5. 


Equations  which  relate  the  left  and  right  eigenvectors  are  derived  from  the  orthogonality 
relationships.  The  first  involves  the  characteristic  matrix  A  and  the  eigenvalues  A. 


It  is  A  U  =  A  (2.42a) 

which  yields  '  ^  u  '  '"I  ^ 
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Making  use  of  Eqs.  (2.41a)and  (2.41b)  in  (2.43), 

AVJU^-*-  AVJ‘U^  =  A  (2.43a) 

so  it  must  be  that  ^u  ^  ®  (2.43b) 

T 

The  orthogonality  relationships  V  U  =  1  can  be  written  in  several  forms.  Since  the  U 
and  V  ntatrices  are  rectangular  of  size  n  x  2n,  the  submatrix  products  of  interest  may  be  n  x  /t  or 
2nx2n. 


The  forms  which  result  in  2n  x  2n  product  matrices  are 

U  =I  and  u'*'v  =  I  (2.44) 

and  the  submatrix  product  matrices  are 

vJ'Uy  +vj‘u^  =  l  (2nj(2n)  (2.45a) 

and  Uy  V  ^  =  1  (2,i  x  2«)  (2.45b) 

The  forms  which  result  in  n  x  n  product  matrices  are 


In  detail. 


U  v'^  = 


and 


vu’^= 


uv'f= 

I  or  vu’’ 

=  I 

u  v'f 

u  u 

I 

0 

U,  vT 

/  u 

v,yj 

0 

I 

V 

u  u 

V  vj 

U  / 

I 

0 

v.u’' 

1  u 

V 

0 

1 

(2.46) 

(2.47) 


(2.48) 


If  the  eigenvectors  occur  in  complex  conjugate  pairs,  there  is  sometimes  an  advantage  to 

iding  the  upper  an( 
complex  conjugate  of  U 


subdividing  the  upper  and  lower  matrices  accordingly  with  an  overbar,  as  for  ,  meaning  the 


up 


u  = 

’  u 

u 

‘  u  u 

up  up 

.  . 

— 1 

s 
_ 1 

(2.49) 


The  four  submatrices  are  each  //  x  n  and  can  be  inverted  because  the  eigenvectors  are  assumed 
here  to  be  independent.  The  diagonal  eigenvalue  matrix  A  must  be  rearranged  accordingly. 


A  = 


A  0 
P 

0  A 

P 


(2.50) 
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2.2.4  Conversiun  from  state  space  to  quadratic  form  and  implications. 

Using  the  upper  and  lower  eigenvector  matrices,  it  is  possible  to  return  to  the  quadratic 
form  of  the  characteristic  equations.  The  derivations  are  done  without  inverting  the  diagonal 
matrix  A,  which  may  be  singular. 

For  the  right  eigenvectors,  using  Eq.  (2.39b),  is  eliminated  from  Eq.  (2.39a)  to  give 
the  quadratic  form. 

lU^A^  +  Cj  U^A  +  KjU,  =  0  (2.51) 

Now,  post-multiply  (2.5 1 )  by  A  and,  since  A,  a  quadratic  equation  in  terms  of  is 

lU  A^  +  C,  U  A  +  K,U  =0  (2.52) 

u  1  u  1  u 


Doing  the  corresponding  process  for  the  left  vectors,  first  pre-multiply  Eq.  (2.41a)  by  A, 

T 

and  then  use  Eq.  (2.4 lb)  to  eliminate  , 

I  A^  +  A  V'^'C,  +  K,  =  0  (2.53) 

u  u  i  u  1 

Taking  the  transpose  yields 

IV  A^  +  CT  V  A  +  kT  V  =0.  (2.53a) 

u  1  u  I  u 

Comparing  Eqs.  (2.52)  and  (2.54),  it  is  seen  that  the  are  involved  with  the  eigenvalue 

T  T 

problem  of  the  transpose  matrices  and  Kj .  If  the  Cj  and  Kj  are  symmetric,  then  V  ^ , 
and  satisfy  the  same  quadratic  equation.  Thus  the  corresponding  eigenvectors  may  be  equal, 

or  differ  at  most  by  a  scalar  multiplier,  which  can  be  different  for  each  mode.  An  equation 
involving  only  is  not  as  simple,  as  shown  in  Eq.  (2.53b). 

-  K’’^  A^  +  cIk'^Vi  a  +  =  0  (2.53b) 

If  the  rectangular  eigenvector  matrices  are  subdivided  into  two  square  submatrices,  which 
are  complex  conjugates  of  each  other,  Eq.  (2.51)  becomes 


+c,ru. 

![_  Ip 


A  0 
P 

0 


P  -I 


•‘■K  ®] 


and  the  results  are  the  n  simultaneous  equation  which  state  the  eigenvalue  problem  of  the  second 
order  differential  equations 


A^  +C,U,  A  +K,  U,  =  0 

P  1  /p  P  1  'P 


(2.56a) 
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and  +  5,^=0  (2.56b) 

The  Eq.  (2.56b)  is  just  the  complex  conjugate  of  Eq.  (2.S6a).  since  the  plant  matrices  are  real. 

The  orthogonality  equation  in  Eq.  (2.57)  cannot  be  further  simpliHed. 

VU’^=rv  U'^  +  V  U'*'1=I  (2.57) 

u  u  |_  up  up  up  upj  ^  ’ 

Example  2.1 

This  example  is  for  a  4x4  state  matrix  with  viscous  damping.  The  purpose  is  to  she  ^he 
numerical  form  of  the  individual  matrices  A  Ap  and  the  product  matrices  in  the 

spending  quadratic  equation. 


I  0 

0  0 

0.4  -0.2  5 

-4 

0  2 

0  0 

♦ 

-0.2  0.8  -4 

4 

M  = 

K  = 

0  0  - 

5  4 

5-4  0 

0 

0  0 

4  -4^ 

-4  4  0 

0 

[-0.4  0.2 

-5  4 

(-0.1376  +  0.52891)  0 

0 

0 

p 

1 

o 

2  2 

0  (-0.2624  +  2.5745  i) 

0 

0 

: 

A  = 

1  0 

0  0 

0  0 

(-0.1376-0.5289  i) 

0 

0  1 

0  0 

0  0 
u 

0  (-0.2624-2.5745  i)  J 

r  (-0.1376+  0.5289  i) 

0 

A  = 

P 

0 

(-0.2624  +2.5745  i) 

_  r  0.5144  +  0.0047  i  -0.8572  +0.0139  i 
"  0.6064  -  0.0020  i  0.3646  +0.0163  i 


The  individual  matrices  in  the  quadratic  are  2x2  and  the  numerical  values  are  listed 

below. 


U,  K-  +C,  U,  A  +K,  U.  =0. 

/p  p  1  /p  p  1  tp 


U,  A 
Ip 


-0.1334  -0.0761  i  5.6414  +1.0670  i 
-0.1584  -0.0877  i  -2.3697  -0.5997  i 
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C,  U,  A 

1  /p 


-0.0128  40.0444  i  0.1032-1.07111 
-0.0256  40. 1013  i  -0.0740  40.5948  i 


K,  U.  = 

l  /p 


0. 1463  40.03181  -5.7446  40.0041 1 
0. 1 840  -0.0 1 35  1  2.4437  40.0049 1 


2.2^  Special  cases  of  combinations  of  and  Kj . 

Special  Case  1:  :  Cj  and  are  symmetric. 

In  general,  if  Cj  =  M'*  C  and  Kj  =  M**  K  ,  then  neither  Cj  or  Kj  will  be  symmetric, 
even  though  C  and  K  are  symmetric.  An  exception  would  be  when  M  is  equal  to  a  scalar  a 
times  the  unit  matrix,  that  is  M=  a  I . 


Another  possibility  for  assuring  that  Cj  and  are  symmetric  is  to  make  a  change  in 
coordinates,  by  weighting  the  original  coordinates  by  the  square  root  of  the  mass  matrix.  The 
result  is  that 

=  andKj=M‘‘^KM‘‘^  (2.58) 

T  T 

which  will  be  symmetric  if  C  and  K  are  symmetric.  For  this  special  case,  Cj  =  Cj  and  Kj  = 
Kp  so  it  is  seen  by  comparing  Eqs.  (2.53)  and  (2.51)  that  the  and  satisfy  the  same 
equations.  Thus,  it  could  be  that  « but  it  is  more  likely  that  they  would  differ  by  a  scalar 

factor  because  of  the  normalization.  The  following  example  demonstrates  this  effect 

Example  2.2 

This  example  shows  that  the  eigenvalues  in  and  differ  only  by  a  scalar  multiplier 
when  Cj  =  Kj  =  0. 


2 

0 

0 

0 

0.5  -0.2 

6 

-3 

0 

2 

0 

0 

♦ 

-0.2  0.4 

-3 

7 

K  = 

0 

0 

-6 

3 

6  -3 

0 

0 

0 

0 

3 

-7 

-3  7 

0 

0 
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-0.25 

0.1 

-3 

1.5 

(-0.0652  +  1.3138  i)  0  0 

0 

0.1  - 

0.2 

1.5 

3.5 

A  = 

0 

(-0.1598  +  2.1777  0  0 

0 

1 

0 

0 

0 

0 

0  (-0.0652-  1.3138  0 

0 

0 

1 

0 

0 

0 

0  0  (-0.1598-2.1777  0  J 

u 

p 


^up 

L  ^/pJ 


-0.0245  +0.7091  i 
-0.0424  +0.6003  i 

0.5394  -0.0081  i 
0.4574  +0.0096  i 


-0.1078  +  0.9954  i 
0.0567  -  1.1758  i 

0.4583  +0.0159  i 
-0.5389  +0.0135  i 


V 

P 


'  -0.0062  -0.4108  i  0.0073  -  0.2 103  i 
0.0073  -0.3484  i  0.0062  +  0.2473  i 

0.5378  -0.0492  i  0.4580  -  0.0278  i 
0.4593+0.00371  -0.5391  +0.04451 


The  result  Is  that  but  the  individual  vectors  are  related  by  scalar  factors.  Adding 

another  subscript  for  eigenvector  number,  V^pj=  (-0.7616  1)  and  =  (-0.4589 1) 

Note  that  two  of  the  four  vectors  are  given,  to  save  space.  The  other  two  vectors  are  the  complex 


conjugates  of  the  two  that  are  given. 


Special  Case  2:  :  =  0  and  is  unsymmetric. 

For  this  special  case,  the  control  forces  are  zero  and  there  is  no  damping.  It  would  not  be 
unusual  for  to  be  unsymmetric.  Then,  from  Eqs.  (2.39), 

.KjU^=U^A  (2.59a) 

and  Ujj  =  A .  (2.59b) 

The  quadratic  equation  reduces  to 

-  KjU^  =  U,A^  (2.59c) 

From  Eqs.  (2.41)  obtained  from  the  characteristic  equation 

V^  =  AVj  (2.60a) 

and  -Vjj  Kj  =  A  vj.  (2.60b) 

It  is  seen  in  Eq.  (2.6()a)  that  the  lower  left  vector  is  related  simply  to  ,  which  is  not  the 
case  when  Cj^O . 
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One  quadratic  equation  is 

-V^K,=A^V'*' 
u  1  u 

(2.60c) 

or 

and  that  with  the 

.k|v,=  v,a^ 

(2.60d) 

The  form  of  X. ,  X.=  -  C  (i).  ±  i  co.  a  /  1-  so  C  =  0  and  X. ,  X.=  ±  i  o). .  The  eigenvalues  of 

j  J  J  M  J  J  J  J 

T 

K|  are  the  same  as  those  for  K^.  but  the  left-hand  vectors  are  still  different  from  the  right-hand 
vectors. 

Special  Case  3:  C^  =  0  and  is  symmetric. 

If  Kj  is  symmetric  with  no  damping,  the  upper  portions  of  U  and  V  satisfy  the  same 

quadratic  equations  as  the  lower  portion  of  U,  but  the  lower  portion  of  V  does  not  This  is  seen 
by  inspection  of  Eqs.  (2.51)  through  (2.54).  This  means  that  the  eigenvectors  are  "equal'*  but 
they  may  differ  by  a  scalar,  depending  on  how  they  are  normalized  or  where  **  m-  =  /  r  ^  r  ■ 

u„  (2.61) 

=U,  A  (2.62) 

V,A.V^_a2  =  U^  =  II,  A  (2.63) 

The  form  of  the  orthogonality  relations  U  =  1  in  terms  of  submatrices  is 

V'^U  +V7u,  =  l  (2.64) 

a  a  l  I  '  ' 
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Example  2  J 

In  this  example,  Cj  =  0  and  Kj  is  symmetric.  If  K  is  synunetric,  then  M  must  be  of  a 
special  form  such  as  M  =  P I ,  where  P  is  a  scalar,  for  to  be  symmetric. 


The  eigenvalues  are 


X,=  1.3150  i,  =  =  \  and  \  =  X2- 

The  eigenvectors  are  in  the  same  order  as  the  eigenvalues. 

‘  0.7095  i  0.9983  i  -0.7095  i  -0.9983  i  ' 
0.6010  i  -  1.1785  i  -0.6010  i  1.1785  i 
0.5395  0.4571  0.5395  0.4571 

0.4571  -0.5395  0.4571  -0.5395 


V 


-0.4103  i  -0.2093  i  0.4103i  0.2093  i  ' 
-0.3476  i  0.2470  i  0.3476  i  -0.2470  i 
0.5395  0.4571  0.5395  0.4571 

0.4571  -0.5395  0.4571  -0.5395  . 
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2.2.6  Recovery  of  the  C I  and  matrices  from  the  eigenproperties. 

The  A  matrix  can  be  recovered  from  the  eigenproperties  by  pre-multiplying  Eq.  (2.42a) 

T 

by  U  and  post-multiplying  by  V  ,  to  produce  Eq.  (2.65), 

UAV'‘’=A  (2.65) 

which,  written  out  in  terms  of  the  submatrices  is 


U  AV'^  U  AVjl  [-C,  -K, 

u  u  u  /  _  1  1 

U.AV''’  U.AVT  I  0 

/  u  /  /  J  L 


(2.66) 


2.2.7  Nonsymmetric  control  matrices,  forced  vibration  problem. 

To  do  a  forced  vibration  problem  with  symmetric,  or  unsymmetric  control  matrices,  G 
and  H,  the  matrix  equations  are 


’m  0  1  fz 

< 

(t)  f  [C+Gl 

-  + 

(K+Hjl  |i(t)l  Jl 

"(t) 

(2.67) 

.0  -IK+HlJ  [z 

(t)  L  [K+HI 

oj  U(t)J  1 

0 

The  advantage  of  using  this  special  form  of  the  redundant  equations  is  that,  if  the  control 
matrices  are  symmetric,  the  closed  loop  M  and  K  matrices  are  again  symmetric.  The 
shorthand  fom.  of  Eq.  (2.67)  is 

M*  T)  +K*  ^  =  F*(t)  (2.68) 

or  q-A*q=[M*r‘  F*(t)  (2.69) 

The  equations  are  decoupled  by  substituting 

n  (t)  =  U  q(t)  (2.70) 

into  Eq.  (2.69)  and  pre-multiplying  by  V  . 
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3.  SENSITIVITY  OF  EIOENPROPERTIES  TO  INDIVIDUAL  ELEMENTS  IN  THE 
DAMPING  AND  CONTROL  MATRICES. 

An  essential  part  of  control  of  vibration  of  structures  is  the  study  of  the  sensitivity  of  the 
eigenproperties,  namely  damping  ratios,  natural  frequencies,  and  mode  shapes,  to  elements  in  the 

damping  and  control  matrices.  The  differential  equations  areMz  +Cz  +  Kz  =  -Gz-Hz, 
where  G  and  H  are  control  matrices.  The  viscous  damping  matrix  C  is  achieved  through 
passive  mechanical,  electromechanical,  or  other  devices  whose  forces  depend  on  the  relative 
velocities  between  two  points  on  a  structure  or  between  one  point  and  ground.  It  is  assumed  that 
there  is  more  freedom  in  designing  the  active  control  forces  on  the  right  side  of  the  equation,  so 

that  it  is  possible  to  apply  a  force  G.^  z.  at  /  proportional  to  the  velocity  at  J ,  for  example.  This 

chapter  deals  with  the  sensitivity  of  eigenvalues  and  eigenvectors  to  changes  in  individual 
elements,  primarily  in  the  C  or  G  matrices. 

An  alternate  form  of  the  differential  equation  is  M  z  +  [C+G]  z  +  [K+H]  z  =  0,  from 
which  the  closed  loop  eigenvalues  are  determined.  So,  an  approach  to  optimum  control  is  this: 
achieve  as  much  damping  as  possible  and  practical  through  optimizing  the  parameters  in  the  C 
matrix  and  then  add  what  more  is  needed  in  the  G  matrix.  The  idea  is  to  design  the  C+G  matrix 
and  the  K+H  matrix,  rather  than  the  individual  matrices,  keeping  in  mind  how  they  are 
physically  achieved. 

The  publications  that  are  closely  related  to  the  present  study  are  summarized.  The 
optimization  of  viscous  damping  to  minimize  free  vibrations  has  been  studied  by  Neubert  [N4], 
Gilheany  [Gl],  Silverberg  [SI],  Mcloughlin  [Ml],  and  Inman  and  Andry  [S],  who  give 
additional  references.  Neubert  [N6j  also  presented  an  approach  for  optimum  sizing  and  location 
of  viscous  damping  to  minimize  forced,  random  vibrations.  Methods  of  active  control  are 
summarized  in  the  book  by  Meirovitch  [M4].  The  need  for  passive  damping  in  feedback 
controlled  flexible  structures  was  emphasized  by  von  Flotow  and  Vos  [FIJ.  Venkayya  and 
Tischler  [VI]  discussed  the  effect  of  frequency  control  on  the  dynamic  response  of  structures. 
Sensitivity  analysis  as  related  to  modal  control  is  the  subject  of  the  book  by  Porter  and  Crossley 
[PI]  in  1972,  which  summarized  the  methods  that  were  developed  to  that  time  and  presented 
some  practical  examples  for  aircraft  stabilization  systems,  cascaded  vehicles,  and  economic  and 
manufacturing  systems.  Somewhat  unique  are  the  sections  dealing  with  systems  which  have 
confluent,  or  repeating,  eigenvalues  rather  than  distinct  eigenvalue  with  linearly  independent 
eigenvectors,  as  assumed  herein. 
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The  contributions  of  the  present  study  that  are  summarized  in  this  chapter  are: 

1 .  The  current,  classical  methods  for  finding  eigenvalue  and  eigenvector  sensitivities  are 

summarized. 

2.  A  direct  mathematical  relationship  is  presented  between  the  coefficients  in  the  complete 

modal  series  method  and  Nelson’s  individual  mode  method  for  Hnding  sensitivities  of 
eigenvectors. 

3.  A  new  method  for  finding  eigenvalue  sensitivities  is  presented,  based  on  derivatives  of 

the  invariants  associated  with  the  similarity  transformations  involved  in  the  eigenvalue 
solutions. 

4.  It  is  demonstrated  that  increasing  damping  may  sometimes  be  deleterious,  because  some 

modes  may  usurp,  or  hog,  the  damping  with  the  result  that  damping  ratios  actually 
decrease  in  other  modes. 

5.  Several  examples  are  given,  that  give  an  engineer  a  practical  feeling  for  the  sensitivity 

of  damping  ratios  and  natural  frequencies  to  changes  in  plant  and  control  matrices. 


3.1  Eigenvalue  sensitivity. 

In  this  section,  methods  are  summarized  in  detail  for  finding  the  sensitivity  of  the 
eigenvalues.  It  is  assumed  that  the  system  is  underdamped  and  the  eigenvalues  occur  in  complex 

conjugate  pairs.  A  typical  pair  of  eigenvalues  is  X,  ^  ^  “r  ^  *  ^Dr  '  l^e 

derivative  of  this  pair  of  eigenvalues  with  respect  to  input  parameters,  the  derivatives  of  the 
modal  damping  ratios  and  natural  frequencies  o)^  are  readily  determined. 

The  eigenvalue  derivatives  are  determined  in  two  ways.  In  the  first  method,  which  was 
developed  during  the  present  study,  the  eigenvalue  derivatives  are  formulated  from  the 
invariants  in  the  similarity  transformation.  Since  the  input  parameters  and  the  eigenvalues  appear 
explicitly  ,  these  derivatives  aid  greatly  in  understanding  the  interaction.  The  second,  standard 
way,  which  is  more  convenient  for  numerical  computations,  is  to  take  the  derivative  of  the 
characteristic  equation.  The  derivative  then  involves  the  eigenvectors  and  the  derivative  of  the 
characteristic  matrix. 

The  matrix  equations  are 

Mz+Ci+Kz  =  -Gi-Hz  (3.1) 

where  the  n  x  n  matrices  M,  C  and  K  are  symmetric  but  G  and  H  may  be  unsymmetric.  The  Eq. 
(3.1)  is  multiplied  by  the  inverse  of  M  to  produce 

z  +  C,  z  +  K|  z  =  -  M'‘  [G  z+  H  z]  (3.2) 

where  Cj=M''c  and  Kj=M'’k  .  The  matrices  Cj  and  Kj  are  usually  unsymmetric.  The 
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resulting  closed  loop  equations  are  written  in  state  space  form  as 


Iy(0-Ay(t)  =0 


(3.3) 


where 


A  = 


-M'‘lC+Gl 

I 


-m'ik+hj 

0 


The  solutions  of  the  homogeneous  equations  are  of  the  form 

y(t)  =  u . 

and  the  eigenvalue  problem  to  be  solved  is 

[-XI  +  a]u  =0. 


(3.4) 

(3.5) 

(3.6) 


From  Eq.  (3.6)  the  eigenvalues  are  determined  and  are  here  assumed  to  be  distinct,  with  no 
repeated  eigenvalues.  The  right  eigenvectors  u  ^  satisfy  the  relationship 

[-X^I  +  a]u^  =0.  (3.7) 

Since  A  is  unsymmetric,  there  is  also  a  set  of  left  eigenvectors  associated  with  the  eigenvector 
problem  involving  A  ,  which  is  the  transpose  of  A.  The  eigenvalues  are  the  same  as  in  Eq. 


(3.7),  but  the  left  eigenvectors  v  ^  satisfy 

[-X^I  +  a'*']v^  =0,  (3.8) 

which  can  also  be  written  in  transpose  form 

v"^^  [•  \  I  +  a]  =  0 .  (3.9) 

The  orthogonality  relationships  between  left  and  right  eigenvectors  are 

v^u^  =0  forr^s  (3.10a) 

and  V  ^  =  1  for  r  =  s  (3. 10b) 

or,  combining  (3. 10a)  and  (3.  lOb),  V^U  =  1  (3. 10c) 

where  U  and  V  are  2n  x  2n  matrices  of  right  and  left  eigenvectors.  It  can  be  shown  that 

V’^AU  =  A  (3.10d) 

where  A  is  a  diagonal  matrix  of  the  eigenvalues. 


Here  Eq.  (3. 10b)  is  the  normalizing  condition  for  the  v  ^  vectors,  where  the  u  ^  have  first 
been  normalized  according  to 

M*  u  ^  =  b^ .  (3.11) 
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3.1.1  Sensitivities  uf  the  eigenvalues  using  the  similarity  transformation  invariants. 

Some  of  the  effects  of  individual  elements  in  the  control  matrix  can  be  anticipated  from 
knowledge  of  the  relationship  between  the  matrix  A  and  the  diagonal  eigenvalue  matrix  A  .  The 
matrices  A  and  A  are  similar,  in  that  they  are  related  by  a  similarity  transformation.  The 
invariants  are  those  associated  with  similarity  transformations,  meaning  that  the  matrices  have 
the  same  eigenvalues,  their  traces  are  equal,  and  their  determinants  have  the  same  value.  The 
equality  of  the  traces  relates  the  system  damping  parameters  to  the  modal  damping  ratios  and 
natural  frequencies,  as  follows. 

Tr  (A)=Tr  ( -(C+G  ))  =Tr  (A)=  j  =  I  -2  (3. 12) 

The  fact  that  the  determinants  are  invariant  means  that 

|a|  =  |m-'  K|  =  a)5  0)^  ...  0)2  ,  (3.13) 

2 

so  that,  even  though  the  O)  ’s  change  due  to  changes  in  C+G,  their  product  does  not  change  if 
the  mass  and  stiffness  matrices  are  not  altered. 

Fom  Eqs.  (3.12)  and  (3.13),  relationships  for  the  sensitivities  can  also  be  determined  by 
taking  derivatives  of  both  sides  of  the  equations  with  respect  to  A.^  .  Thus  the  following 

Important  results  are  obtained. 


a(Tr  (A))  _  -a(Tr(C+G))  ^  a(X  ( -2  co^) 

aA..  aA..  aA..  aA.. 

>j  u  y  u 

a|M'‘  k|  a(  0)2  0)^ . . .  0)^ ) 

aA..  aA.. 


(3.14) 

(3.15) 


y  y 

The  Eq.  (3.12)  is  written  in  detail  as  Eq.  (3.16)  and  taking  the  partial  of  (3.16)  with 
respect  to  A.,  confirms  Eq.  (3.33),  which  is  derived  separately  below. 


2n  2n 

l\  -  1 

2n2n 

(3.16a) 

r=l  i=l 

and 

La''' 

r=l  “ 

(3.16b) 
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Example  3.1:  System  with  One  Degree  of  Freedom 


For  the  simple  one-degree-of-freedom  system  with  a  mass  m  supported  by  a  spring  of 
stiffness  k,  a  dashpot  with  damping  rate  c,  and  a  control  force  of  -g  x  the  A  matrix  is 


A  = 


m'Vc+g)  -m‘*k 


1 


and  the  diagonal  eigenvalue  matrix  is 
A  “dI 


0 


0 


(3.17) 


(3.18) 


0  -  CjtOj  -  i  tOpj 

The  following  relationships  are  obtained  by  setting  the  traces  and  determinants  equal. 


-  m''  (c+g)  =  -  2  C,  Wj 


and 


(3.19) 

m  k  =  o)j  +  co^j  =  (Oj .  (3.20) 

Taking  the  partial  derivatives  of  Eqs.  (3.19)  and  (3.20)  with  respect  to  g,  the  results  are 


-I 


2  2  2  2 

.  0).  +  CO...  =  0).  . 


-1 


and 


0  =co. 


dcOj 


(3.21) 


(3.22) 


From  Eq.  (3.22),  if  (0j?t  0,  then  9(Oj/3g  =0  which  shows  that  changing  g  does  not  result  in  a 
change  in  (0^. 

Example  3.2:  System  with  two  degrees  of  freedom 

If  the  system  has  two  degrees  of  freedom,  with  the  damping  matrix  C=0,  and  the  mass 
matrix  is  diagonal,  the  A  and  A  matrices  would  be  4x4,  as  follows. 


A  = 


-G  -K  ■ 

I  0  J" 


-Gjj/mj  -G|2/m|  -  Kjj/mj  - 

-G2j/m2  '^2’^^2 


0 

1 


0 

0 


0 

0 


(3.23) 


and 


A  = 


0 

0 

0 


0 

0 

0 


0 

0 

0 


0 

0 

0 


"2  -I 


(3.24) 
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From  the  traces  and  determinants,  the  following  relationships  arise. 

(K,,K22-K,2K2,)/(mjm2)  =  X,  X,  (a]  (3.25) 

-G|  i/m|-G22/ni2=  *  -  ^1 '  ^2  ‘  ^2“ 

The  fact  that  the  eigenvalues  must  be  real  or  occur  in  complex  conjugate  pairs  is  con^rmed  by 
Eq.  (3.26),  since  the  left  side  is  real.  Further  it  shows  that  if  G|  j=  ^22~  ^  ^12 

system  will  be  unstable  since  and  must  then  be  real  and  of  opposite  signs,  with  and  0)2 

real  and  positive. 


By  taking  the  partial  derivative  of  Eqs.  (3.25)  and  (3.26)  with  respect  to  G^^  ,  equations 
for  the  sensitivities  are  obtained. 


2  at02  2  ^“1 


and 


l/(2m  )=C 


do) 


'8G..  ^2  00 


8o), 


+0) 


•dG. 


+  Q), 


2dG. 


(3.27) 


(3.28) 


II  II  — II  — 11 

An  interesting  .situation  ari.ses  according  to  Eq.  (3.27),  namely  that  a  change  in  G^  ^  could 

produce  changes  in  (Oj  and  (O2  and  if  one  of  them  increases,  the  other  must  decrease. 


3.1.2  Sensitivities  of  the  eigenvalues  in  terms  of  the  eigc.ivectors. 


To  find  the  sensitivity  of  the  eigenvalues  with  respect  to  a  parameter  in  location  A.^. 

take  the  partial  derivative  of  Eq.  (3.7)  with  respect  to  A.j .  The  well  known  equation  [R1,W1] 

for  solving  for  this  derivative  in  terms  of  the  eigenvectors  and  the  derivative  of  A  is  given  in  Eq. 
(3.29). 


dX  dA 

- =  v*  -  u  . 

d  A-j  dA.. 

«J  U 


(3.29) 


The  right-hand  side  of  Eq.  (3.29)  is  a  scalar.  Taking  the  partial  of  A  with  respect  to  A.^ 

results  in  a  matrix  having  unity  in  the  ij  position  and  zeroes  for  all  the  other  matrix  elements.  , 
using  v.^  to  designate  the  \th  element  of  the  vector  v  ^  an  ^  the  j.'/t  element  of  u^,  the  Eq.  (3.29) 

becomes 

dX 

=v.  u.  (3.30) 

dA..  jr 
U 

The  2n  x  2n  sensitivity  matrix  is  now  formed  [PI],  which  shows  how  each  element  A  ^  affects 
the  eigenvalue  X  ^ , 
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r  dx  1  _ 

5  =  ^  =v 


(r.  i,j=  1,2,  ..2n) 


(3.31) 


The  sums  of  sensitivity  matrices  are  of  interest  and,  from  Eqs.  (3.31)  and  (3.10c),  are 

2n 

^S^=VU^=I  (3.32) 

r=l 

where  the  V  and  V  are  2/t  x  2n  matrices  made  up  of  the  right  and  left  eigenvectors  respectively. 

Writing  the  Eq.  (3.32)  explicitly  for  the  sum  of  the  diagonal  elements  of  the  sensitivity 
matrices  shows  a  significant  result  in  Eq.  (3.33),  which  is  the  same  as  Eq.  (3.16b). 


(3.33) 


2n 

y  s..  ^ 

=  y 

Li  iir 

Lu  9  A.. 

r=l 

\  H 

r=l 

2n 

ys..= 

y  = 

Li  Ur 

Li  9  A.. 

r=l 

r=l  “ 

(3.34) 


The  product  of  two  sensitivity  matrices  for  two  different  eigenvalues  is  zero,  as  can  be 
seen  from  the  following,  which  follow  from  Eqs.  (3. 10). 


SS=vu^v  u^  =  0  (r?ts) 

rsrrss  ^ 

SS=v  u^v  u^=S  (r=s) 
rr  rrrr  r 


(3.35a) 

(3.35b) 


If  the  system  is  underdamped,  the  eigenvalues  occur  in  complex  conjugate  pairs,  and  the 
following  relationships  apply; 


+  icOp,  (3.36a) 

X  X  =  0)^  and  X  +  X  =-  2C  o)  .  (3.36b) 

rrr  rr^rr  ^  ' 

If  there  are  2n  eigenvalues,  there  are  only  n  different  modal  frequencies  O)^ .  The  is  the 

modal  frequency  of  damped  free  vibrations. 

By  differentiating  the  two  expressions  in  Eq.  (3.36b)  with  respect  to  A.^ .  the  derivatives 
of  the  modal  damping  ratio  and  natural  frequency,  that  is  0A.^  and  9(0^  9A.^,  are  readily 

determined  in  terms  of  the  derivatives  of  X  and  X  . 

r  r 
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EXAMPLE  3.3:  Ten-bar  truss  with  eight  degrees  of  freedom 


The  ten-bar  Uuss  studied,  and  shown  in  Fig.  3.  1,  is  one  that  has  been  suggested  by  some 
investigators  as  one  of  a  group  of  "standard”  trusses  and  was  used  as  an  example  in  references 
[N4]  and  [N6]  .  It  is  not  necessarily  a  practical  truss,  but  is  used  for  comparison  purposes.  The 
truss  lies  in  the  x-y  plane  and  has  four  moving  nodal  points,  numbered  1,  2,  3,  and  4,  which  are 
free  to  move  in  the  x-  and  y-directions.  There  are  two  anchored  nodal  points,  numbered  S  and  6. 
Each  truss  bay  is  square,  360  inches  on  each  side.  The  cross-sectional  areas  are  given  in  Table 

3.1.  The  material  of  each  bar  is  aluminum  with  a  Young's  Modulus  E=:10xl0^psi.  Half  the  mass 
of  each  bar  was  lumped  at  each  end  of  the  bar  in  the  present  analysis. 


Fig.  3.1  Ten-bar  truss 


Table  3.1.  Areas  of  Truss  Bars 


Bar  No. 

2 

Area(in  ) 

Bar  No. 

Area(in^) 

m 

A 

m 

m 

A 

m 

1 

31.5 

6 

0.5 

2 

0.1 

7 

7.5 

3 

23.0 

8 

20.5 

4 

15.5 

9 

21.0 

5 

0.1 

10 

O.I 

Optimization  of  passive  viscous  damping,  representable  by  dashpots  connected  between 
nodal  points,  was  presented  in  (N4J  and  (N6J  for  this  truss.  In  practice,  the  control  forces 
represented  by  the  Gx  would  be  superimposed  on  the  damping  forces  generated  by  Cx.  In  order 

to  isolate  the  effects  of  the  control  forces,  no  passive  damping  is  included  in  these  first  examples, 
so  the  matrix  C=0  initially. 

It  is  assumed  that  a  control  force  can  be  located  at  any  of  the  four  moving  truss  nodes  and 
that  it  is  proportional  to  velocity  at  that  or  any  other  node,  in  the  x-  or  y-direction.  The  positions 
in  the  G  matrix  follow  the  numbering  of  the  nodes  on  the  truss.  Thus  rows  I,  2  and  3  correspond 
to  lx,  ly  and  2x  on  the  truss,  where  2x  means  the  displacement  at  node  2  in  the  x-direction. 
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Hence  G23“^Iy:2x  ^57“^3x5x  example.  The  sensitivities  of  the  eigenvalues  to 

individual  elements  in  the  G  matrix  were  found  using  Eq.  (3.12).  The  relationships  in  Eqs. 
(3.22)-(3.23)  help  to  understand  and  confirm  the  results,  but  they  were  not  used  directly  in  the 
computations  of  sensitivities. 

The  modal  damping  ratios  ^  of  the  clos^  loop  system  are  of  particular  interest  For  this 

truss,  it  was  found  that  the  danq)ing  ratios  in  different  modes  are  affected  by  various  choices  of 
location  and  magnitude  of  G..  In  Fig.  3.2,  the  damping  ratios  Cg.  C5  and  are  plotted  for  both 

positive  and  negative  values  of  G^^  .  It  is  seen  that  for  small,  positive  values  of  G^^  the 

damping  ratios  and  are  both  positive  and  increase  almost  linearly  with  G^^  .  All  the  modal 

damping  ratios  are  positive  in  this  range  and  the  system  is  stable.  For  larger  values  of  G^^, 

increases  at  a  faster  rate  and  ^5  decreases,  but  they  remain  positive.  The  curves  are 

antisymmetrical  with  respect  to  the  origin. 


Fig.  3.2  Variation  of  anting  with  G^^ 
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Fig.  3.3  versus  for  r  =  6, 5  and  8 

The  sensitivities,  or  derivatives  of  and  with  respect  to  arc  shown  in  Fig.  3.3. 

These  are  related  directly  to  the  slopes  of  the  curves  in  Fig.  3.2.  When  the  sensitivity  is  positive, 
the  damping  ratio  is  increasing.  Thus  for  the  slope  is  seen  to  be  always  positive,  while  for 

the  slope  is  positive  for  small  G^^  and  negative  for  larger  values. 

The  damping  ratios  of  the  same  three  modes  are  shown  in  Fig.  3.4  versus  G^^  .  The 

corresponding  natural  frequencies  are  shown  in  Fig.  3.5  and  the  damped  natural  frequencies 

are  plotted  in  Fig.  3.6  versus  G.^.^.  It  is  seen  that  damping  is  highest  in  the  8th  mode  and  that 

approaches  1.0  or  critical  damping  as  G.^.^  becomes  large,  while  and  as  well  as  the 

damping  ratios  for  the  other  five  modes,  remain  relatively  small.  The  range  of  the  natural 
frequencies  for  small  damping  is  from  131  to  796  rad/s.  When  damped  natural 

frequency  i>ntl  the  modal  frequency  (i)g=  598  rad/s,  so  there  is  a  decrease  of  about  25  %  in 

cOg.  At  the  same  time  the  other  seven  modes  all  show  increases  of  from  0.03  to  7.5  %  for  their 

(I)  's. 
r 
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It  should  be  mentioned  that  it  becomes  difficult  to  trace  the  modes  as  the  modal  damping 
ratio  increases.  The  mode  numbers  are  here  decided  in  association  with  small  control  forces.  If 
the  modes  are  ordered  according  to  the  magnitude  of  (Oj^.  the  damped  natural  frequency,  then 

what  was  the  eighth  mode  for  small  control  forces  becomes  the  first,  or  lowest  frequency,  mode 
for  large  control  forces. 

In  Fig.  3.7,  the  damping  ratios  are  plotted  for  an  input  control  force  for  .  The  sensi¬ 
tivities  are  not  plotted,  but  are  readily  visualized  as  the  slope  of  the  ^  curves.  The  and 
sensitivities  are  positive  for  positive  but  ^  shows  a  negative,  decreasing  trend  in  region. 
Thus  the  system  is  unstable  if  only  the  G.^^  control  force  is  applied.  If  it  is  desired  to  increase  the 
damping  in  modes  5  and  8  and  decrease  that  in  mode  6,  then  a  positive  G^^  should  be  used.  A 
negative  G.^^  would  produce  the  reverse  effect.  The  effect  of  G^^  on  the  modal  damping  ratios  is 
similar  to  that  of  . 

The  relationships  between  and  G..  or  Gjj  are  nonlinear.  The  plots  show  a  nearly  linear 

range  only  for  small  control  forces.  These  nonlinearities  are  due  to  changes  in  mode  shapes. 
Large  viscous  damping  or  velocity-proportional  control  forces  tend  to  stiffen  the  structure  and 
change  the  mode  shapes.  The  eigenvectors  are,  in  general,  complex  and  their  sensitivities  are 
readily  calculated,  as  discussed  in  Section  3.2. 

Because  of  the  nonlinear  behavior,  the  question  arises  as  to  how  well  the  effects  of  G^^, 

G^y ,  and  G^.^  superimpose.  This  is  answered  graphically  in  the  plots  in  Fig.  3.8,  where 

is  shown  for  the  four  control  forces  applied  simultaneously,  with  G^^and  G^^  positive  but  with 

G^^  and  G^^  negative,  to  increase  For  the  dashed  line,  the  separate  effects  were  added,  while 

for  the  solid  line,  the  forces  were  applied  simultaneously.  It  is  clear  that  superposition  applies 
for  small  control  forces  only.  Note  that  the  range  of  the  control  force  on  Fig.  3.8  is  much 
smaller  than  on  the  previous  figures.  As  the  magnitude  of  the  four  control  forces,  applied 
simultaneously,  approaches  3(K)()  lb  s/in,  the  value  of  ^^=1.77  and  the  damping  ratios  of  all  the 

other  modes  are  less  than  0. 1 .  so  mode  6  usurps  the  dissipation. 
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Fig.  3.8  Cg  with  effects  of  Gjj,  and  q^lied  siinultaneously  (solid) 

and  supeiiinposed  (dashed) 

Example  3.4  Effect  of  initial  passive  damping,  ten-bar  truss 

Passive  viscous  damping,  which  is  repiesratable  mathematically  by  dashpots  connected 
in  parallel  to  each  truss  member,  tends  to  produce  dissipation  in  every  mode  of  vibration  and  the 
system  is  stable.  If  there  is  passive  damping  and  only  one  active  control  force  is  added  in  an 
off-diagonal  position,  such  as  G^g,  then  the  resulting  motion  can  be  stable  or  unstable, 

depending  on  the  relative  magnitudes  of  the  passive  and  active  forces.  Using  the  ten-bar  truss  as 
an  example,  proportional  damping  was  included  by  having  a  dashpot  parallel  to  each  truss 

member,  so  the  C  =  P  K,  with  p  =  3. 1623x10*^  s.  Then  the  value  of  G^^g  was  varied  from  -200  to 

200  Ib^  s.  The  result  is  shown  in  Fig.  3.9  as  a  plot  of  ^  and  versus  G^g.  The  damping 

ratios  for  the  other  six  modes  are  not  small,  but  they  are  affected  little  by  G^g.  The  solid  line  is 

for  the  system  with  both  viscous  dancing  and  the  single  control  force.  The  dot-dash  line 
represents  the  values  when  the  viscous  damping  is  zero  and  only  the  single  control  force  is 
active.  For  the  latter  situation,  the  trace  of  C+G  is  zero,  so  ^  and  ttmd  to  have  opposite  signs 

and  the  system  is  unstable.  When  the  passive  damping  is  added,  the  curves  shift  upward  and  the 
system  becomes  more  stable. 
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r=4 


undamped 


T - 1 - 1 - }-— I - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

■  IM  200 

Gjg  (lb  s/in) 

Fig.  3.9  C2  and  versus  Gjg  with  and  without  proportional  passive  damping 


3.1.3  Summary  and  condusions 


1.  The  modal  damping  ratios,  ^  ,  are  odd  functions  of  the  G-  ,  so  their  sensitivities 
d^dG.^  are  even  functions.  If  there  is  no  passive  damping,  the  ^  tend  to  be  positive  when 
the  diagonal  elements  G..  are  positive.  However,  the  ^  can  be  positive  or  negative  when 
the  off-diagonal  elements  G.^  are  positive. 

2.  The  C.  are  nonlinear  functions  of  the  various  G..  .  For  small  control  forces,  the 

’r  ij 

relationship  is  nearly  linear.  The  nonlinearity  for  larger  G.^  is  due  to  a  changes  in  mode 
shapes. 

3.  In  view  of  the  changes  in  mode  shapes,  complex  eigenvectors  must  be  used  and 
proportional  damping  cannot  be  assumed  except  for  small 

4.  The  study  emphasized  the  effects  of  individual  elements  G.^  in  the  G  matrix.  If 

combinations  of  the  G.^  are  used,  corresponding  to  mutiple  control  forces,  their  individual 

effects  superimpose  accurately  only  for  small  control  forces.  This  is  because  the  changes  in 
mode  shapes  for  individual  G..  are  not  the  same  as  for  combinations  of  G.. . 

5.  As  has  been  well  known  previously,  the  addition  of  passive  damping  tends  to  have  a 
stabilizing  effect  in  combination  with  the  types  of  control  forces  considered  herein. 
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3.2  Sensitivity  of  the  eigenvectors. . 

The  sensitivity  of  eigenvectors  to  changes  in  the  plant  matrices  is  of  interest  when  modal 
solutions  are  used.  Some  of  the  instances  where  they  are  important  are  as  follows: 

1.  In  predicting  changes  in  kinetic  or  potential  energies,  the  modal  derivatives  appear  ex¬ 
plicitly. 

2.  When  a  perturbation  technique  is  used  for  pole  placement,  it  may  be  advantageous  to 
predict  the  change  in  eigenvectors. 

3.  If  proportional  damping  is  assumed,  then  the  corresponding  assumption  is  that  the  mode 
shapes  do  not  change.  The  eigenvectors  derivatives  indicate  how  much  change,  or  error,  is 
acceptable  in  the  plant  matrix  parameters  in  order  to  be  able  to  assume  that  changes  in  mode 
shape  are  negligible. 


An  important  point  is,  even  if  mode  shapes  representing  displacements  do  not  change,  the 
eigenvectors  in  State  Space  do  change  if  the  eigenvalues  change.  This  is  because  the  upper  right 
vectors,  for  example,  are  equal  to  the  lower  vectors  (the  mode  shapes)  times  the  eigenvalue 
matrix,  even  for  non-symmetric  plant  matrices. 

3.2.1  Sensitivity  of  the  eigenvectors  by  the  complete  modal  series  method. 


du. 

The  eigenvector  derivative  is  written  in  short-hand  as  u..^  .  For  the  derivative  of 


the  ith  vector,  start  by  taking  the  derivative  of  [A  -  X.  I]  u.  =  0 . 

[A  -  li  u.,  +  [A,^. ,  -  X. ,  .  ,1)  u.  =  O 

I  *  j  ak/  ak/  I  ak/  ‘  i 

It  is  known  that  the  derivative  of  the  eigenvalue  is  found  from 


\  ’ak/  '"i 


Sk/“i 


(3.37) 


(3.38) 


It  appears  that  the  eigenvector  derivative  might  be  found  by  simply  inverting  the  [A  -  X.  I]  matrix 

in  Eq.  (3.37),  but  that  matrix  is  singular  and  cannot  be  inverted,  so  another  approach  must  be 
found.  Two  possibilities  are  summarized.  First,  a  common  approach  when  working  with 
eigenvectors  is  to  represent  a  solution  in  terms  of  a  weighted  sum  of  the  eigenvectors  [PI],  so 
that  u.,  is  expressed  as 


u.. 


2n 


,  =y  Y  u  =UY-- 

ak/  ■  , '  «j  j  ii 
j=l 


(3.39) 
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The  Eq.  (3.39)  requires  that  aii  the  be  known.  Substituting  Eq.  (3.39)  into  (3.37)  gives 


J=* 


(3.40) 


To  solve  for  the  y  j. ,  pre-multiply  Eq.  (3.40)  by  v^  and  make  use  of  the  relationships 


v^A  u.  =0  n^j 

n  j  J 

v^A  u.  =  X.  n=j 

n  j  j  J 

Thus,  only  the  jr/j  term  in  the  series  is  non- zero,  with  the  result 


(3.41a) 

(3.41b) 


(X.  -  X.)  y  .  =  v^  A,^. ,  u. 
J  t  'ij  n  ak/  1 


(3.42) 


V.  A,„,,  u. 

y .  =  -  ^  — -  fori^tj. 

9  (X.  -  X.)  •' 


(3.43) 


For  i  =  j,  the  denominator  of  Eq.  (3.43)  goes  to  zero,  so  y..  must  be  found  from  a  different 
equation. 

It  is  interesting  to  repeat  the  derivation  more  efficiently  by  substituting  the  matrix 
representation  of  Uj,  of  Eq.  (3.39)  into  Eq.  (3.37)  and  pre-multiply  by  V^,  the  full  matrix  of 
left  eigenvectors. 


vV-A,Il  UYi^-vTA.^  «, 


V'^A  U  =  A  =  ‘  ^ 


and  suppose,  for  example,  that  i=2.  Then  Eq.  (3.44)  becomes 


'>12  A.^ 


■  1^21 


and  Eq.  (3.43)  follows. 


To  determine  y..  we  may  use  the  normalizing  condition 


(3.44) 


(3.45) 


(3.46) 


lx  2n  2ii  X  2ii  2m  x  I 

M  u.  =  b. 

1  I  I 


and  take  its  derivative  with  respect  to  a^^^ 


(3.47) 
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This  gives  M  Uj  +  u[  ilj  +  u>  =  0 

If  M  is  symmetric,  the  equation  can  be  written 

.  u.  +  2  uT  M  u.  . ,  =  0 
I  ak/  I  1  i.ak/ 

and  the  derivative  of  the  vector  is  replaced  by  its  series  form 


or 

yielding,  finally,  for 


•*«  ^11 

u.  u.  +  2u.  MY  Y  U  =0 
1  akf  1  1  jti  'J  J 

M,  ,  u.  +  2  b.  Y-  +  2  M  y  y  u.  =  0 
I  ak/  I  I  *11  I  •  , '  ij  j 

J=} 


Having  tnc  y.  ,  the  eigenvector  derivative  may  be  calculated  from  Eq.  (3.39). 


(3.48) 

(3.49) 

(3.50) 

(3.51) 

(3.52) 


Before  leaving  this  derivation,  it  is  worthwhile  to  contemplate  the  expression  in  Eq. 
T 

(3.47),  namely,  u.  M  u.  =  b,  .  One  simplification  would  occur  if  the  vectors  were  normalized 

T 

so  that  bj  =1.  A  benefit  would  also  accrue  from  prudent  choice  of  M,  such  that  M  u.  =0. 

For  symmetric  plant  matrices,  this  would  occur  if  IV1=  M 

3.2.2  Eigenvalue  sensitivity  by  Nelson’s  single  mode  approach. 


The  disadvantages  of  the  complete  series  approach  are  that  ail  the  modal  vectors  must  be 
known  and  summing  them  for  each  eigenvector  derivative  is  tedious. 

In  developing  Nelson’s  |N2|  single  mode  approach,  substitute  Eq.(3.38)  into  (3.37)  to 

obtain 


I''  •  \  "  “i’  ak/  =  •''■ak/  "i  "i  1  “i  (3.53) 

=  F 

I 

Rather  than  use  a  modal  series  to  represent  the  solution.  Nelson  represented  u.,  ^  as  the 
sum  of  a  homogeneous  and  a  particular  solution,  namely 

“i’ak/ 

The  contribution  of  Nelson  is  that  the  vector  h.  may  be  calculated  directly,  without  knowing  the 
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remaining  eigenvectors.  Suppose  that  u.,  ^  is  represented  by 


“raw  =  « 


+  a..  <  u. 


(3.55) 


where  the  sub  1  and  2  just  mean  portions  of  the  vectors  h.  and  u. .  Then,  the  special  form  of 
the  vector  h.  is  solved  for  by  zeroing  out  the  row  and  column  of  the  [A  -  X.  I]  matrix.  Since 
the  rank  of  the  reduced  [A  -  X.  I|  must  be  2n-l,  the  remaining  (2n-i)  x  (2n-l)  matrix  may  be 
inverted.  The  appearance  of  the  partitioned  matrices  is  as  follows. 


tA-X.I|,j  0 


10  |[|j 

0  IA-X,IJ„  I": 


=  0 

F, 


(3.56) 


L  "31  i‘^33J  ‘  '• 

which  is  easy  to  program  on  the  computer.  By  inserting  the  digit  1  in  location  jj  on  the  main 
diagonal  of  the  matrix  on  the  left,  the  matrix  is  restored  to  size  2n  x  2n  and  h^  =  0  is  assured  by 

zeroing  the  F  vector  in  the  same  row.  The  choice  of  which  vector  to  designate  as  the  }th  row  is 
arbitrary,  but  to  avoid  numerical  difficulties  Nelson  suggested  finding  the  row  where  the  compo¬ 
nents  of  the  vectors  Uj  and  v.  are  large  by  investigating  the  products  of  the  absolute  values 

abs(u.)  abs(v..). 

The  a.,  is  found  from  the  normalizing  condition,  Eq.(3.1 1),  after  taking  the  derivative  as 


in  Eq.  (3.49)  above 


M,,^. ,  u.  +  2  uT  M  u. 


u^M,  u.  +  2ul*^Mlh. +  a..u.  1  =0 
so  a.  =-  ~  I  uJ  M,  . ,  u.  +  2  u7  M  h.  I 

II  2b.  \  ‘  ‘  ‘  / 

Here  a  simplification  will  result  if  IVI  is  chosen  such  that  =0. 


(3.49) 

(3.57) 

(3.58) 
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3.2J  Reiationships  between  the  complete  series  approach  and  the  single  mode  method. 

In  using  the  complete  series  representation,  from  Eq.  (3.39),  the  eigenvector  derivative  is 
given  by 

J=l 

By  Nelson’s  single  mode  approach,  the  derivative  is  given  by  Eq,  (3.45) 


u.  „i.i  =h. +  a..u. 
I  ak/  I  11  I 


(3.45) 


It  is  possible  to  represent  h.  as  a  series  summation  of  all  the  modal  vectors,  weighted  by  a  factor 


lij  •  “ 


It  follows  that  the  series  coefficients  are  related  by 


(3.59) 


Y.  =Ti..  fori^j 

'ij  'ij  J 

Y.  =Ti.. +  a..  fori=j. 

'll  'll  II  ■' 


(3.60a) 

(3.60b) 
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Example  3^.  Eigenvector  sensitivity  by  the  complete  series  approach 


The  matrices  in  this  example  are  given  in  state  space  form. 


M*  = 


’  1 

0 

0 

0  ‘ 

’0 

0 

5 

-4 

0 

2 

0 

0 

K*  = 

0 

0 

-4 

4 

0 

0 

-5 

4 

5 

-4 

0 

0 

0 

0 

4 

-4 

-4 

4 

0 

0 

so 


A=  [M*]-'  K*  = 


0-5  4 

0  2-2 
0  0  0 
1  0  0 


(3.61) 


(3.62) 


TT  ♦ 

The  right  eigenvectors  were  normalized  according  to  U  M*  U  =I .  The  sensitivities  calculated 
are  those  of  the  right  eigenvectors  with  respect  to  A,-  or  u,  .  The  right  eigenvectors  arc: 


U  = 


-0.3645  -0.3645  -0.6059  -0.6059 

-0.4284  -0.4284  0.2577  0.2577 

0.6672  i  -0.6672  i  0.2341  i  -0.2341  i 

_  0.7843  i  -0.7843  i  -0.0996  i  0.0996  i  J 

The  eigenvalues  are  Xj.  ^2^  0.5463  i ;  X^=  ±  2.5887  i . 

The  sensitivities  "■2’ai2 


(3.63) 


Xj,ai2=<J1562  X2,^, 2=0.1562  X^,^^2=-0A562 


The  values  of  are  listed  in  matrix  form,  their  subscripts  determining  their  position  in 


the  matrix. 


0  0.1429i  -0. 1 27 li  0.0828 i 

-0.1429  i  0  -0.0828  i  0.1271  i 

-0.0460  i  -0.0300  i  0  -0.0302  i 

0.0300  i  0.0460  i  0.0302  i  0 


(3.64) 


The  partial  derivatives  of  the  right  vectors  are  arranged  in  columns 


u 


i,al2 


du. 


d\ 


12 


-0.0253  i 

-0.0727  i 

0.1445 

0.0912 


0.0253  i 

0.0727  i 

0.0145 

0.0912 


-0.0460  i 

-0.0248  i 
0.0036 
0.0156 


0.0460  i 
0.0248  i 
0.0036 
0.0156 


(3.65) 
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Example  3.6  Eigenvector  sensitivity  by  the  single  mode  approach. 


The  plant  matrices  and  the  eigenvectors  are  the  same  as  in  Example  1.  The  matrices  for 
finding  the 


*^1  = 


■'^’aI2“! 


+  lv 


^’al2  “i  ^  “l 


0.3715 
-0.0669  I 
0.1042  i  I 
0.1225  ij 


The  [A  -  Xj  I]  matrix,  after  modification  of  the  third  row  and  third  column  is 


(3.66) 


Modified  |A  -  X.  I|  = 


-0.5463  i 


0 


0 

0 

0 


Modified  Fj  = 


-0.5463  i 
0 

1 

0.3715 
-0.0669 
0 

0.1225  i 


0 

0 

1 

0 


4 

-2 

0 

0 


(3.67) 


(3.68) 


Now  the  matrix  h.  is  found  by  multiplying  the  inverse  of  the  modified  lA-XjI]  matrix  times  the 


modified  Fj  matrix. 


Thus 


Since  andbj=1.0 


0.1042  i 

0.1654  i 
0 

0.0786 


and 


ttjj  =  -  u|m  hj=  0.2166  i 

0.0253  i 
0.0727  i 
-0.1445 
-0.0912 


“i,aI2  =  »'r 


(3.69) 


(3.70) 

(3.71) 


which  agrees  with  the  first  column  in  the  Eq.  (3.65)  above,  obtained  by  the  complete  series 
method. 


-38- 


4.  OPTIMIZATION  OF  PASSIVE  DAMPING 

4.1  Optimization  of  viscous  damping. 

4.1.1  Achieving  uniform  modal  decay  rates  for  free  vibration  problems. 

In  the  present  context,  passive  viscous  damping  is  that  representable  by  dashpots,  which 
may  be  connected  to  ground  or  between  two  points  on  a  structure. 

Mz+Cz+Kz=0  (4.1) 

In  free  vibration  problems,  it  is  sometimes  desirable  for  the  motion  of  each  mode  to  decay  at  the 
same  rate.  This  may  be  achieved  with  mass-proportional  damping,  such  that  C  =  2  a  M  .  The 
displacements  are  represented  in  terms  of  modal  coordinates. 

z(t)  =<l>q(t)  (4.2) 

T  T 

Then,  if  O  M  0=  I ,  the  product  d>  C  <1>  =a  I  and  a  typical  equation  for  the  ith  mode 
is 


The  eigenvalues  are 
so  the  decay  rate  is  the  same  in  each  mode. 


q  +2  a  q  +  (0^  q 
=  -  a  ±  i  0)^  -a^ 


=  0. 


=  -  a  ±  i  0), 


(4.3) 

(4.4) 


If  the  mass  matrix  is  diagonal,  of  size  n  x  /i,  then  n  dashpots  are  theoretically  required.  If 
the  system  is  three-dimensional,  a  dashpot  is  required  at  each  mass  in  each  direction,  and  one  end 
must  be  connected  to  an  effective  ground.  Thus  achieving  mass-proportional  damping  may  not 
always  be  practical. 

For  stiffness-proportional  damping,  C=  P  K  and  d>  =  p  [to^]^  and  the  decoupled 
equations  are  of  the  form 

2  •  2 

q  +Po}  q+0)q  =0. 

^  r  T  r  ^  r 

The  standard  form  of  the  decoupled  equations  in  terms  of  the  modal  damping  ratio  ^  is 

q^  q^  +  (dJ  q^  =0.  (4.5) 

so  the  equivalent  modal  damping  ratio  is  found  from 

2 ;  (0  =P  0)^ 

which  shows  that  the  higher  modes  decay  fastest,  unless  the  ^  ^  are  designed  so  they  decrease 

with  frequency.  Practically,  stiffness-proportional  damping  may  be  easier  to  achieve  than 
mass-proportional  damping  because  mechanical  devices  can  be  connected  between  nodes,  rather 
than  to  ground.  But  then  active  control  may  need  to  be  added  to  provide  sufficient  dissipation  in 
the  lower  frequency  modes. 
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4.1.2  Critical  damping  fur  matrix  equations 


The  idea  of  a  critical  damping  matrix  was  presented  by  Inman  and  Andry  [12]  in  1980. 
Usually  the  mass  matrix  M  and  the  stiffness  matrix  K  are  positive  definite,  although  K  may  be 
positive  semi-definite  for  a  free  structure  which  possesses  rigid  body  modes. 

Now  change  variables  by  letting  x  =  z. 

X  X  +M'‘^KM‘'^x=  0 

X  +Cj  X  +Kj  X  =  0 

Now  critical  damping  is  defined  in  reference  (I3I  as  Ci=  2  kJ^  . 

The  following  cases  are  then  possible: 

1 .  If  C I  =  1 1  the  system  is  critically  damped  and  each  eigenvalue  is 
a  repeated  negative  real  number.  All  the  eigenvectors  are  real. 

2.  IfCj-C|  is  positive  definite  then  the  system  is  overdamped  and  each 
eigenvalue  is  a  negative  real  number.  All  the  eigenvectors  are  real. 

3.  IfCj-C|  is  negative  definite  then  the  system  is  underdamped  and  the 
eigenvalues  occur  in  complex  conjugate  pairs,  each  with  a  negative  real  part 

4.  lfC|-t!|  is  indefinite  or  sign  variable,  then  mixed  damping  occurs.  At 
least  one  mode  will  be  overdamped  and  at  least  one  mode  will  be  underdamped. 

Related  references  and  concepts  were  summarized  in  the  paper  by  Liang,  Tong,  and  Lee 
[Lljwho  give  definitions  for  strongly  complex  and  weakly  complex  eigenvectors.  Caughey  and 
O’Kelly  1C2]  gave  a  necessary  and  sufficient  condition  for  decoupling  of  the  Eqs.  (4.6). 


(4.6a) 

(4.6b) 
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4.13  Optimizatiun  of  modal  damping  ratios,  ten-bar  truss  example. 

4.13.1  Formulation  of  the  optimization  problem. 

In  a  practical  truss,  bars  are  connected  at  pinned  joints.  The  viscous  damping  force  on 

pin  j  may  be  F.  =  (  z.  -  z  )  if  the  dashpot  is  between  stations  j  and  j+1  .  If  the  dashpot  is 

con  nected  between  station  J  and  ground,  then  F^  =  z. .  It  is  expensive  to  connect  damping 
devices  between  each  pair  of  nodes.  So  the  questions  to  be  addressed  are: 


1 .  Which  modal  damping  ratios  are  to  be  maximized? 

2.  If  a  limited  number  of  dashpots  are  to  be  used,  where  should  they  be  connected 
on  the  structure? 

3.  Once  the  best  locations  are  determined,  how  much  damping  should  be  used  at 
each  location? 

For  the  underdamped  situation,  the  complex  eigen  pairs  are  written  as  in  Eq.  (4.7),  with  ^ 
modal  damping  ratio. 

X.,X|  =  -C.<o,±ja.,^l-  Cf  (4.7) 

The  optimization  problem  is  specified  by  defining  appropriate  objective  functions. 
Suppose  it  is  desired  to  maximize  the  damping  ratios  of  certain  modes.  Then  the  objective 
function,  OBJ,  could  be  specified  as 

OBJ  =  J]  Oj  Cj  (4.8) 

I 

where  the  are  arbitrary  weighting  factors  to  be  chosen  as  positive  real  numbers  or  zero.  In 
addition  a  bound  is  put  on  the  total  amount  of  damping  to  be  used.  This  is  done  by  defining  a 
constraint  function  G(  I ),  using  weighting  factors  and  as  the  viscous  dashpot  constants. 

G(l)  =  y  Y  C  -C_  =  0.0  (4.9) 

'  'in  m  T 

m 
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4.U.2  Natural  frequencies  and  mode  shapes  for  proportional  damping. 


The  ten-bar  truss  of  interest  lies  in  the  x-y  plane  as  shown  in  Fig.  3.1,  page  23  .  The 
connections  are  assumed  to  be  frictionless  pins.  The  nodes  S  and  6  are  anchored  and  the  other 
four  nodes  are  free  to  move,  with  x  and  y  components  of  motion.  The  bar  numbers  are  in 


brackets.  The  bar  material  is  aluminum  with  =  10  x  10^  psi  for  each  bar.  The  dashpots  are 
not  shown,  but  each  dashpot  is  connected  between  the  same  two  nodes  as  a  parallel  bar,  with  the 
same  element  number,  so  is  connected  parallel  to  bar  m.  The  axial  stiffness  of  each  bar  is 


A  E 

— .  The  bar  cross-sectional  areas  are  listed  in  Table  4.1.  The  truss  was  originally 


m 

designed  to  carry  heavy  vertical  loads  at  nodes  2  and  4  .  Consequently,  the  bars  2,  S,  6  and  10 
have  a  much  smaller  cross-sectional  area  than  the  oth^r  six  bars.  The  truss  bays  are  square,  360 
inches  on  each  side. 


The  stiffness-proportional  dashpot  constant  is  arbitrarily  defined  as 


C  = 

m 


Pm 


^  ■ 

where  is  an  factor  to  be  chosen.  If  P^,  is  a  constant,  then  the  damping  is  stiffness 

proportional.  For  the  first  example,  all  P^  =  10,  The  mass  matrix  is  diagonal,  with  half  the 
mass  of  each  bar  lumped  at  the  node  at  each  bar  end.  The  damping  ratios  and  natural 
frequencies  for  this  situation  are  given  in  Table  4.2.  The  damping  factors  for  the  higher  modes 
are  the  largest.  The  natural  frequencies  range  from  131.13  to  796.21  rad/s  and  are  well 
separated. 


The  mode  shapes  for  zero  damping,  or  for  proportional  damping,  are  shown  in  Figs.  4.1. 
For  modes  1  and  3,  the  shapes  are  somewhat  like  lateral  bending  modes  of  the  entire  truss.  For 
the  higher  frequency  modes,  the  stretching  of  the  lighter  bars,  namely  bars  2,  6  and  10, 
dominates  the  motion. 


Table  4.1  Member  areas,  lengths,  stiffnesses,  and  damping  factors  for  P^=  1.0  . 


m 

A  (in') 
m 

L„(in) 

K^^(lb/in) 

C^^(lb-s/in) 

1 

31.5 

360,000 

875  (KK) 

276.699  30 

2 

0.1 

360.01M) 

2  778 

0.878  41 

3 

23.0 

360.(K)0 

638  889 

202.034  41 

4 

15.5 

360.(MK) 

430  556 

136.153  62 

5 

0.1 

360.IMM) 

2  778 

0.878  41 

6 

0.5 

360.(KK) 

13  890 

4.392  05 

7 

7.5 

509.117 

147  314 

46.584  75 

8 

20.5 

509.117 

402  658 

127.331  65 

9 

21.0 

509.117 

412  479 

130.437  30 

10 

0.1 

509.117 

1  964 

0.621  13 
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Tabic  4.2  Damped  modal  damping  ratios  and  natural  frequencies  for  i.o. 


i 

(0 .  (rad/s) 

1 

0.0207  3350 

131.1301 

2 

0.0274  1571 

173.3921 

3 

0.0426  4408 

269.7048 

4 

0.0515  5221 

326.0448 

5 

0.0730  5049 

462.0119 

6 

0.0948  6770 

599.9960 

7 

0.1051  3561 

664.9360 

8 

0.1258  9253 

7%.2143 

4. 1.3.3  Sensitivities  uf  modal  damping  ratios  and  natural  frequencies. 

The  first  derivatives,  dCldC  were  calculated  by  the  methods  of  section  3.1.2,  after 

I  iti 

determining  the  derivatives  of  each  eigenvalue  d  X.  /  dC^  and  its  complex  conjugate  d  Xj  /  dC^. 
The  derivatives  of  the  damping  ratios  are  tabulated  in  Table  4.3.  There  are  80  entries  in  the 
table,  since  there  are  8  modes  and  10  bars.  The  total  of  each  row  and  column  is  also  given  and 
the  sub-totals  add  to  an  overall  total  of  10.27867.  There  is  a  big  variation  in  the  magnitudes  of 
the  sensitivities,  which  is  the  reason  for  listing  values  to  six  decimal  places.  The  list  of  dashpot 
locations  in  order  of  sensitivities  is  2,  10, 6,  5, 4, 9,  3, 8,  1, 7.  By  bar  stiffness,  from  smallest  to 
the  largest,  the  order  is  10,  2,  5,  6,  7,  8,  9,  4,  3,  1  .  Thus,  there  is  a  direct  correlation  between 
sensitivities  and  stiffnesses,  with  largest  sensitivities  going  with  dashpots  parallel  to  bars  with 
lowest  stiffnesses.  Location  7  seems  to  be  most  out  of  order,  but  notice  that  there  are  two  nega¬ 
tive  values  in  the  column  for  bar  7.  In  a  ranking  of  potential  effectiveness  according  to  the 
absolute  sum  of  the  columns,  the  value  for  da.shpot  7  would  be  0.037637  and  it  would  precede 
dashpot  3  in  the  sensitivity  list. 
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MODE  2. 


MODE  I .  •*»/* 


MODE  3.  •*»/* 


MODE  4.  m.«4  lAO/s 


MODE  5,  **o/s 


MODE  6,  *00.00  «*o/t 


MODE  7.  ••••«  ■«/* 


MODE  8.  7t*.],  iAO/S 


Fig.  4. 1  Mode  shapes  of  the  undamped  ten-bar  truss. 


Table  4.3  Sensitivities 


X  10‘  at  p  =  1.0  . 
dC  •" 

m 


i 

m=l 

m=2 

m=3 

m=4 

m=5 

m=6 

1 

0.002  025 

0.004  389 

0.000  900 

0.001  138 

0.007  674 

0.000  035 

2 

0.000  197 

0.021  051 

0.001  518 

0.000001 

0.154  036 

0.000  238 

3 

0.000  180 

0.004  650 

0.005  255 

0.010  230 

0.006  373 

0.000  355 

4 

0.000  054 

4.336  IW 

0.000  037 

0.000  021 

0.000  486 

0.000  238 

5 

0.008  992 

0.045  043 

0.004  783 

0.001  471 

0.007  432 

0.001  360 

6 

0.011  616 

0.(K)9  531 

0.005  680 

0.000  325 

0.000  890 

0.069406 

7 

0.000  203 

0.018  862 

O.(KX)  531 

0.000  367 

0.000  079 

2.122  320 

8 

0.007  902 

O.IHK)  946 

0.017  544 

0.051  749 

0.000  016 

0.007  902 

0.023  581 

4.440  572 

0.036  248 

0.065  302 

0.176  986 

2.234  647 

i 

m=7 

m=8 

m=9 

m=10 

All  m 

1 

-0.001  810 

0.(X)6  954 

0.(X)1  747 

0.007  211 

0.030  263 

2 

-0.012  466 

().0(X)  762 

O.(XX)  059 

0.029  303 

0.194  699 

3 

0.006  633 

0.(X)7  399 

0.(X)1  391 

0.045  729 

0.045  729 

4 

0.000  104 

0.0(X)  117 

0.000  000 

1.853  090 

6.223040 

5 

0.003  230 

0.014  265 

0.012  294 

0.039  581 

0.138  451 

6 

0.003  364 

0.003  603 

0.031  427 

0.081  899 

0.217  741 

7 

O.OCC  344 

O.IXK)  039 

0.(X)1  667 

1.164520 

3.308  932 

8 

0.009  686 

O.IXX)  024 

0.010  797 

0.020  834 

0.119  812 

0.009  085 

0.033  163 

0.061  254 

3.197  829 

10.278  67 

4.U.4  Modal  usurping  uf  damping;  variation  of  C3  and  with  Cj . 

As  the  values  of  are  increased  and  the  damping  becomes  nonproportional,  the  mode 

shapes  change  and  the  eigenvectors  become  complex  numbers,  so  the  eigenvalue  derivatives 
change.  The  variation  of  and  with  is  shown  in  Fig.  4.2.  All  the  other  were  taken 

as  zero.  The  plot  shows  that  for  0  <  Cj  <  3,  both  C5  and  increa.se  almost  monotonically. 
However,  for  Cj>  3,  the  rale  of  increa.se  of  becomes  larger  while  that  of  decreases.  The 
slopes  of  these  lines,  which  are  the  sensitivities,  are  plotted  in  Fig.  4.3,  where  it  is  seen  that 
/9Cj  changes  sign  and  finally  approaches  zero  as  Cj  increases.  So  mode  6  usurps  the  increased 
damping  and  mode  5  damping  approaches  zero,  becau.se  of  changes  in  mode  shapes. 
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and 


0  2  *  6  a  lO  I*  ■«»  "o 

Fig.  4.2  Variation  of  Cj  and  ^  with  Cj. 


Fig.  4.3  aCj/SCj  and  . 
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4.U.5  Optimization  of  C3  with  respect  to  and  . 

Now  we  return  to  the  optimization  problem  stated  at  the  beginning  of  this  chapter.  The 
goal  is  to  optimize  the  damping  ratio  of  mode  3  by  placing  dashpots  parallel  to  bars  3  and  4, 
which  are  horizontal  members  at  the  bottom  of  the  truss.  The  table  of  sensitivities  shows  us  that 
the  C3  is  quite  sensitive  to  dashpots  placed  in  these  positions. 

The  objective  function  is  taken  as  C3  with  the  design  variables  and  in  the 
constraint  function  G(l). 

OBJ  =  C3 

G(1)  =  C3-»-C4-  0.0 

An  advantage  to  limiting  the  design  variables  to  two  is  that  a  two-dimensional  plot  can  be  shown 
of  the  interaction  between  C3  and  the  two  design  variables.  Contours  for  ^3  =  0.05,  0.10,  0.15 

amd  0.20  are  shown  in  Fig.  4.4.  The  straight  lines  are  constraint  lines  which  were  chosen  to  be 
approximately  tangent  to  the  given  contour  lines.  Thus  each  constraint  line  corresponds  to  a 

different  value  of  Cj.  in  G(l).  The  slopes  would  be  -45*^  if  the  vertical  and  horizontal  scales  were 

equal.  Since  the  contour  lines  are  convex  toward  the  feasible  region,  the  optimum  solution, 
which  maximizes  ^3,  will  be  along  the  constraint  boundary  line. 
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The  results  are  summarized  in  Table  4.4  for  three  values  of  C.^. ,  namely  Cj,  =  1157.3, 
1838.5  and  2525.4  Ib-s/in.  The  associated  values  of  optimum  are  close  to  0.10,  0.15  and 
0.20.  It  is  interesting  to  note  that,  even  though  the  increments  in  are  constant  at  0.05,  the 
corresponding  increments  in  and  are  not  constant,  indicating  the  nonlinear  nature  of  the 
optimization  problem. 

Table  4.4  Optimum  values  of  with  and  the  design  variables. 

C4  C3  C.J. 

160.928  996.373  0.0994  46  1157.301 

719.330  1119.184  0.1470  48  1838.514 

1240.957  1284.475  0.2015  59  2525.432 
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4.2  Representation  and  uptimization  of  viscoelastic  damping  in  structures. 

4.2.1  Fitting  curves  to  data  fur  complex  moduli  of  viscoelastic  damping  materials. 

A  great  deal  of  data  exists,  in  chart  form,  for  the  Young’s  Modulus  and  loss  factor  of 
commercially  available  materials  that  are  used  for  damping  structural  vibrations.  Typically  all 
the  data  is  presented  in  one  nomograph  as  a  function  of  frequency,  reduced  frequency  and 
temperatures,  it  is  more  convenient  to  represent  the  data  in  equation  form,  using  only  three  short 
equations  and  seven  parameters  for  each  material. 

The  basic  equations  are  not  new,  but  a  method  of  fitting  is  suggested  which  is  better  than 
currently  used  in  the  outside  literature.  A  step-by-step  procedure  is  outlined  for  finding  the 
seven  parameters.  As  an  example,  specific  values  are  given  for  the  parameters  for  Soundcoat 
601,  which  has  a  transition  temperatures  of  75“  F.  In  reference  |N7|,  curve-fitting  is  also  done 
for  Soundcoat  606  and  609,  which  have  transition  temperatures  of  125  and  175®  respectively. 
The  fit  obtained  is  demonstrated  in  plots  of  Young's  modulus  E,  loss  factor  q,  and  loss  modulus 
Et)  which  compare  the  fitted  curves  with  experimental  data  which  was  supplied  by  the 
manufacturer. 


List  of  Symbols 


P 

Eu 

E/ 

E(to,T) 


E*(a},T) 

E(o),T)Ti(a),T)  = 
f 


n 

ri(o),T) 

T 

T 
'  () 

T 


Shift  factor  depending  on  temperature  T 

Used  for  initial  frequency  shift  at  T=To  . 

Constant  in  equation  for  o.^ . 

Maximum  value  of  E,  occurring  at  high  frequencies 
Minimum  value  of  E,  occurring  at  low  frequencies 
Real  part  of  Young's  Modulus 

Complex  modulus  =  E(a),T)  1 1  +  j  T|(to,T)] 

Loss  modulus 

0)  /27t  =  frequency  (Hz). 

Shape  and  magnitude  parameter  for  E(a),T)  and  Ti(a),T) 

Loss  factor 

Variable  temperature 

Transition  temperature  for  material 

Temperature  parameter  in  equation  for  O.^, . 

Frequency  (rad/s). 
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4.2. 1 . 1  Equations  and  parameters 

There  are  three  basic  equations  for  the  curve  fitting,  as  given  in  Ref.  [Nl]  .  They  are 
expressions  for  the  Young's  Modulus  E(ci),T),  the  loss  factor  ilCco.T),  and  the  temperature  shift 
factor  Oj. . 

E(a),T)  =  E#  +  Eu  I  1 - - - K  1 

n  n  Eu(  P(»xx.j.)" 

"  2E(co,T)[l+(Pa)a.j.)"J^ 

C|(T-To) 

(T-TJ 

There  are  seven  parameters  involved,  namely  E^,  E^,  Cj,  T^,  T^,  P  and  n.  It  can  be 
seen  in  Eq.  (4.10)  that  the  vertical  position  of  the  E(a),T)  curve  is  controlled  by  E^  and  E^. 
Further,  when  (i)=(),  E(a),T)  =  E/  and  when  o)  is  large,  E(to,T)  =E^.  The  shape  of  the  curves  and 
also  the  magnitude  of  ri(a),T)  are  directly  controlled  by  n.  The  ccj.  and  p  are  factors  of  frequency 
(I)  and  have  the  effect  of  shifting  the  curves. 

Note  that  when  T  =  Tq  the  value  0^=1,  so  there  is  no  temperature  shift.  The  temperature 
Tq  is  called  the  transition  temperature  of  the  material  and  the  maximum  value  of  T|(ti),T)  usually 
occurs  near  that  temperature.  Thus  the  P  may  be  used  to  locate  the  E(to,T)  in  the  correct  lateral 
position  on  a  frequency  scale,  when  T=To.  Then  at  other  temperatures,  the  tXj  automatically 
accounts  for  the  temperature  changes,  if  it  is  chosen  correctly.  Since  the  Oj  in  Eq.  (4.12) 
depends  greatly  on  two  parameters,  Cj  and  T^,  finding  the  best  combination  is  not  always  easy. . 
The  temperature  T^  is  usually  chosen  somewhere  in  the  glassy  region  and  the  value  of  Cj  often 
lies  between  1  and  12. 

In  the  book  by  Nashif,  etal  INI]  it  is  suggested  that  some  of  the  fitting  can  be  done  by 
knowing  where  the  loss  modulus  E(a),T)ri((0,T)  reaches  its  maximum.  From  Eqs.  (4.10)  and 
(4.1 1),  the  loss  modulus  is 


(4.10) 

(4.11) 

(4.12) 


E(a),T)  Ti(a),T)  = 


n  Jt  Eu(  Ptoot.^)" 

2  [l  +( PoKx/J^ 


Taking  the  derivative  with  respect  to  O)  yields 


(PoKXj)'^  =1  at  Pto  for  max.  loss  modulus. 


(4.13) 


(4.14) 
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Using  this  result  in  Eq.  (4.13),  the  maximum  value  of  the  loss  modulus  is  obtained. 

[E«a,T)n(«.,T)]^,=^  (4.15) 

It  is  suggested  that,  if  the  other  parameters  are  known,  then  the  value  of  n  may  be  calculated 
from  Eq.  (4.  IS).  A  practical  problem  arises  in  using  this  relationship  for  Soundcoat  601  because 
the  maximum  of  the  loss  modulus  occurs  at  a  frequency  higher  than  the  range  covered  by  the 
given  composite  plot.  This  was  overcome  by  extending  the  frequency  scale  to  f=1.0E8  and  by 
extending  the  straight  temperature  lines  to  the  higher  frequency  range. 

It  is  easier  to  find  the  maximum  value  of  the  loss  factor,  which  occurs  near  the  transition 
temperature  and  is  easily  identified,  as  can  be  seen  in  Fig.  4.5  .  Hence,  the  derivative  of  the  loss 
factor  ■ri((D,T)  is  needed  with  respect  to  o).  Taking  the  derivative  and  setting  the  result  equal  to 
zero  gives  the  following  relationship. 

'  u 

It  is  convenient  to  take  the  square  root  and  define  the  resulting  specific  value  of  (jkoOy)”  where 
the  maximum  occurs  as  R. 


(pOJOj)"  = 


E/ 

E/  +Eu 


(4.17) 


The  expression  for  the  maximum  value  of  loss  factor  is,  in  terms  of  R,  and  using  Eqs.  (4.10)  and 
(4.11), 


Ti(w,T,a_) 


n  Ji  Ey  R 


T  max  R 

2jE,  +  E  -  Ti+R1 

[  '  u  I  +R  L  J 

^4 


(4.18) 


For  many  practical  materials  R  is  small  compared  to  1,  or  R  «  1  because  E^  «  E^. 
Taking  advantage  of  these  inequaliiie.s,  Eq.  (4.18)  becomes 


T|(0),T,0Lj.) 


max 


n  7t  Ey  R 


2{e,.E„r} 


(4.19) 


If  it  is  also  true  that  E/ «  E^  R,  then  Eq.  (4.18)  may  be  simplified  to  Eq.  (4.20),  but  Eq.  (4.18)  is 
preferred  for  its  accuracy. 


n  Jr  Eu  R 


n  71 

T" 


(4.20) 
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The  Eqs.  (4.16)  '(4.20)  do  not  appear  in  Ref.  INIJ  and  have  been  developed  by  the 
writer.  It  appears  that  Eq.  (4.20)  may  be  used  easily  to  find  the  value  of  n,  but  since  n  also 
controls  the  "bandwidth"  of  the  TKtoT.Op)  curve,  judgment  is  deferred  until  a  rit  of  actual  data  is 
attempted  below. 

4.2. 1.1  Reading  data  from  a  composite  plot  or  nomograph. 

The  data  for  Soundcoat  601  is  shown  in  Fig.  4.S.  There  are  three  vertical  scales.  On 
the  right  the  value  of  frequency  f  in  Hz  is  the  ordinate;  on  the  left,  the  values  of  E  and  T].  On 
the  bottom  horizontal  scale  the  value  of  oc^f  may  be  read.  This  value  may  be  used  to  find  o.^. 

In  obtaining  data  from  the  plot  for  a  certain  temperature  and  frequency,  the  procedure  is 
to  start  with  a  combination  of  frequency  f  and  temperature  T,  for  example  say  f=1000  Hz  and 

T=75‘’  F.  So  draw  a  horizontal  line  through  f=10^  Hz  and  find  its  intersection  with  the  slanting 

temperature  line  where  T=75‘’  F.  Then  drop  a  vertical  line  from  that  intersection  point  and  mark 
the  crossing  points  on  the  E  and  T)  curves.  The  values  of  E  and  t)  are  then  read  on  the  left  hand 
scales.  But  now  an  important  point  is  this:  if  the  vertical  line  is  projected  down  to  the  horizontal 
scale,  the  value  of  fot.^  may  be  read  and,  knowing  f,  the  value  of  Oj  is  readily  calculated.  For 

this  example  fa.^.  =  10^  Hz,  so  o.^  =  1  because  both  T  and  T^  =  75”  F.  If  T=50”F,  then  the 
value  of  Oj  =27,  and  the  same  value  should  result  no  matter  what  value  of  frequency  f  is  taken 
for  this  T. 

In  Table  4.5,  the  values  of  E,  T)  and  Eti(o),T)  are  listed  for  Soundcoat  601  for  a 
temperature  of  75‘’F. 


Table  4.5  Dynamic  Properties  of  Soundcoat  601 

601 


f(Hz) 

E(psi) 

tlE(psi) 

T  =  75°F 

n 

1 

210 

57 

0.28 

10 

360 

180 

0.50 

100 

1  260 

945 

0.75 

1000 

6  900 

6  900 

1.00 

4000 

18  000 

14  4(K) 

0.80 

1.0E4 

30  9(K) 

18  450 

0.60 

4.0E4 

90  000 

32  4(K) 

0.36 

1.0E5 

180  000 

43  2(K) 

0.24 

4.0E5 

210  (K)0 

37  800 

0.18 

1.0E6 

240  (K)0 

31  2(H) 

0.13 

4.0E6 

270  (K)() 

21  6(K) 

0.08 

1.0E7 

300  000 

19  5(K) 

0.065 

4.()E7 

300  000 

12  OIK) 

0.04 
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Shear  Storage  Modulue 


Teaperature  T,  (*F) 


Reduced  Frec|uency  f0(^ 


T  (®F) 


Fig.  4.5  Nomograph  for  G  and  11  as  functions  of  f,  otj. .  and  f  a^, . 
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4.2.1.3  Step>by*step  procedure  for  fitting  equations  to  data. 

Now  the  following  procedure  for  determining  the  parameters  is  proposed. 

1.  Read  the  values  of  E/  and  E  from  the  curve  of  E  vs  f  •  For  601  the  values  are 
E/=18Upsi  and  Eu  =  300,000  psi. 

2.  Read  the  value  of  curve.  Use  this  to  calculate  n  from  Eq.  (4.20).  For 

601,  =  l.(X)  and  the  value  of  n  =  0.637.  To  check  using  Eq.  (4.18),  R=0.02449 

and  n=0.669,  which  results  from  the  more  accurate  calculation. 

3.  On  the  loss  factor  curve,  at  the  transition  temperature  Tq  ,  read  the  value  of  f  OLp  at 
which  the  loss  factor  is  maximum.  Since  o.^  =  1  at  the  transition  temperature,  the 
value  of  f,n  =  I  (KM)  Hz  for  601. 


or 


4.  Read  the  value  of  E  at  frequency  (%  =  2jtf,n  at  temperature  Tq  .  This  value  of  E 
will  be  used  to  calculate  fi,  the  shift  factor,  when  o.^  =  1. 


E((D,n,T,)  =  E  +  E  I  1 - - - K  I 

I  “  L  I  Kn)  J 

(Po))"  = 


Fu+E  /  -  E((0,„  ,T^) 


(4.21) 

(4.22) 


For  601,  E(a),n  ,T  ^)  =  69(M)  psi,  so  =  0.022913  and  P=5.63E-07. 


5.  Next,  account  for  the  effect  of  the  variable  temperature  T,  by  determining  oCp  Since 

the  data  for  601  is  plotted  versus  f  the  value  of  a.j.  may  be  read  directly  for  each 

frequency  f  at  any  temperature  T.  The  values  so  determined  are  listed  in  Table  4.6  for 
601.  From  any  two  pairs  of  values,  T  and  cCp  the  values  of  Cj  and  T^  may  be 

calculated  for  use  in  Eq.  (4. 12),  which  may  then  be  used  to  find  the  value  of  (Xj.  for  any 

temperature. 

Table  4.6.  Values  of  oup  versus  Temperature  for  Soundcoat  601 


T  ("F) 

«T 

0 

1.0  E7 

25 

3.0  E3 

50 

2.7  El 

75 

l.O 

100 

8.0  E-2 

125 

1.2  E-2 

150 

2.0  E-3 

175 

4.3  E-4 
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4.2.I.4  Results  of  curve-fitting  for  Soundcoat  601. 

The  values  of  Oj  at  the  highest  and  lowest  temperature  listed  in  Table  4.6  were  used 
to  determine  C,  and  T  ,  which  are  listed  in  Table  4.7.  These  values  of  Ci  and  T_  were  then 

I  oo  •  oo 

used  to  calculate  the  o.^.  at  all  temperatures. 

Table  4.7  Listing  of  Parameters  for  Curve-fitting  for  Soundcoat  601 


T 

0 

75°  F 

Eu 

300,000 

E/ 

180 

(Eq)inax 

43,200 

at  f= 

100,000 

'Imax 

1.00 

at  f,n= 

1000 

R 

0.02449 

P 

5.630  E-07 

n 

0.669 

T 

OO 

-1 15.75°  F 

Cl 

9.788 

The  correcmess  of  P  is  best  checked  at  the  transition  temperature,  when  Oj.  =1.0.  In 
Figs.  4.5  and  4.6.  plots  of  E,  t|.  and  En  are  depicted  versus  f  for  601  for  three  temperatures  T=50, 
75  and  125"F.  In  Fig.  4.5,  the  comparison  of  the  manufacturer's  data  (dotted)  and  Eq.  (4.10) 

fit  (solid)  is  especially  good  for  E(a})  at  the  transition  temperature  To=75”F  and  in  the  frequency 
range  near  1000  Hz.  This  shows  that  the  value  of  p  is  well  chosen.  The  shape  of  the  curve  at 
higher  and  lower  frequencies  is  not  as  good  and  depends  in  part  on  choice  of  n.  At  the  other  two 

temperatures,  T=50  and  1 25  ‘’F  the  fit  in  the  center  region  of  the  curves  is  excellent  indicating 
that  the  Oy  is  about  correct  Again  the  fit  is  not  as  good  at  the  ends  of  the  curves. 

Both  the  maximum  value  of  r|  and  the  shape  of  the  r|(a),T)  curve  depend  on  n.  Plots 
of  n  versus  frequency  f  for  the  three  temperatures  are  shown  in  Fig.  4.6.  The  maximum  value 
matches  well  and  the  shape  of  the  curve  is  reasonable,  but  the  fit  near  the  ends  of  the  curves  is 
not  as  good.  On  the  log-log  plot,  the  shape  of  the  q  vs.  f  is  symmetrical  and  covers  several 
decades  in  magnitude  while  the  experimental  data  is  unsymmetrical  and  covers  little  over  one 
decade  in  magnitude. 
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pin  I  I  ~i — r 


Fig.  4.6  E  vs.  f  for  601  at  T  =  50, 75,  and  125**  F,  experimental  and  analytical. 


4.2. 1.5  Impruving  the  fit  by  modificatiun  of  the  equatiun  for  t). 

The  fit  of  11  depends  greatly  on  the  value  of  n,  which  affects  both  the  magnitude  of  ii 
and  the  curve  shape.  A  larger  value  of  n  produces  an  ii  curve  of  sharper  curvature;  a  smaller  n 
flattens  the  ti  curve.  However,  if  the  value  of  n  is  increased  over  those  listed  in  Table  4.7  for  any 
one  of  the  three  example  materials,  then  the  magnitude  of  the  magnitude  of  the  maximum  of  T| 
will  be  too  high.  Hence  it  is  suggested  that  the  Eq.  (4. 1 1 )  for  tKcd.T)  be  modified  to  include 
another  constant  to  adjust  the  magnitude  of  ti(0),T)  as  in  Eq.  (4.23).  Then  the  n  may  be  used 
mainly  to  control  the  curvature  of  effective  bandwidth  of  the  curve.  Unfortunately  the  choice  of 
.1  affects  p. 


Ct,  n,  Jt  Eu(  Po)a^)"« 
2E(a),T)[l  ■^v  pcoa.j.)"»p 


(4.23) 


The  following  procedure  is  therefore  suggested  for  improving  the  fit  of  il(ci),T).  Follow 
the  usual  procedure  as  outlined  in  the  steps  above  to  find  n.  Then  modify  n  slightly  to  improve 
the  fit  and  call  this  value  n,  for  "trial  value  of  n".  Use  this  value  n, ,  in  Eq.  (4.24)  to  find  p.  The 
value  of  is  calculated  simply  by  taking  the  ratio  of  the  n’s  or  =  n  /  n^.  Some  typical 


numbers  for  several  values  of  n,  are  given  in  Reference  IN7J. 


(Pa)OCp)"‘  = 

4.2.i.6  Summary  and  conclusions 


0.5 


(4.24) 


Based  on  the  standard  equations,  a  new  procedure  is  suggested  for  calculating  the 
required  parameters.  The  resulting  parameters  are  listed  in  tabular  form  in  Table  4.7  and  curves 
are  presented  showing  that  the  fits  for  E  and  Er)  are  good.  The  fit  for  T)  is  not  as  good,  but  in 

calculations  using  the  complex  modulus  as  E*  =  E(l+Jr|)  the  value  of  Et|  is  important  and  not 
the  value  of  q  itself.  However,  a  method  is  suggested  for  improving  the  prediction  of  q  and  it 
happens  that  the  prediction  of  E  and  Eq  is  also  improved,  as  demonstrated  in  Figs.  12-14,  over 

the  frequency  range  of  10  to  10^  Hz.  The  predictions  are  made  for  temperatures  from  50  to 

2()0‘’F. 
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4.2  J  Representation  uf  viscoelastic  damping  by  complex  modulus. 

As  has  been  explained,  the  complex  modulus  of  practical  uamping  materials  varies 
with  frequency.  However,  in  the  frequency  range  where  the  loss  factor  is  a  maximum,  the  plot 
of  loss  factor  versus  frequency  is  fairly  flat.  So.  if  the  natural  frequencies  of  the  structure  are  in  a 
narrow  range,  reasonable  estimates  of  the  eigenproperties  may  be  obtained  by  using  the  complex 
modulus  in  that  frequency  range  and  assuming  it  is  constant 

For  simple  structures,  the  element  stiffness  is  proportional  to  Young’s  modulus  of  the 

material.  If  the  Young's  modulus  is  complex,  as  E  ==  E(l+iT))  and  the  natural  frequencies  will 

be  of  the  form  O)^  (  1  +  i  r)^) .  The  eigenvalues  found  in  the  state  space  will  be  of  the  form  of 

A.^=sqrt(-(i)“  (  1  +  i  Ti„))-  'Thu''  the  loss  factor  for  each  mode  is  determined  by  the  complex 

natural  frequencies  and  if  the  structure  has  the  same  loss  factor  throughout,  each  mode  will  have 
the  same  loss  factor  and  the  eigenvectors  will  be  real  »nd  the  same  as  for  the  undamped  structure. 
This  is  demonstrated  in  Example  4. 1 


Example  4.1 


M  = 


I  0 
0  2 


K  =  (l  +0.1  i) 


3(M)  -50 
-50  400 


The  eigenvalues  and  eigenvectors  of 


K  are 


so 


0)^  =  18X.76+i  18.88, 
Tl,=Tl2 


(o]  =311.24  +  i  31.12 

=  0.10 


and 


<t>=\ 


0.3029  0.9530 
0.6739  -0.2142 


The  eigenvalues  and  eigenvectors  in  state  space  are 


A.  =±(-0.6861  +i  13.7562) 
1 


=  ±  (-0.8810  +  i  17.6639  ) 


with 


U. 


0.3029  0.3029  0.9530  0.9530 
0.6739  0.6739  -0.2142  -0.2142 


and  U  =  U,  A 
u  / 


In  Example  4.2,  the  loss  factors  vary  from  0.05  to  0.10  for  the  elements  in  the  stiffness 
matrix.  As  a  result,  the  loss  factors  for  the  two  modes  are  found  to  be  0.054  and  0.096,  as  shown 
in  the  Example.  In  addition,  the  mode  shapes  are  complex. 

Example  4.2 


M  = 


1  0 
0  2 


K  = 


The  eigenvalues  and  eigenvectors  of  M'  K  are 


300(1  +i0.10)  -50(1  +i0.08) 

-50(1  +i0.()8)  40()(  1  +i().()5) 
-I 
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(i)|  =  188.86+110.19,  0)^  =311.14  +  129.81 

so  n,  =0.054,  Tij  =0.0% 

0.2988-10.0263  0.9547  +  10.0082  ' 
0.675 1  + 1 0.0058  -0.2 1 1 3  + 1 0.0 1 86 

From  the  state  space,  the  eigenvalues  are 

=  ±  ( -0.3705  +  1  13.7476)  ^3  4  =  ±  (-0  8441  + 1  17.6594  ) 

with  !?2?2]  ^ 


and  <I>  =[<p  I  2  ]  ~ 
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4.2.3  Overdamped  mini-oscillatur  approach  fur  representation  of  viscoelastic  dements. 

As  mentioned  above,  the  problem  with  mathematical  representation  of  viscoelastic 
damping  is  that  the  material  properties,  particularly  the  elastic  m^uli,  vary  with  frequency. 
Since,  for  free  vibrations,  each  mode  vibrates  at  its  natural  frequency,  a  different  elastic  modulus 
should  be  used  for  each  mode.  But  that  does  not  result  in  a  typical  eigenvalue  problem.  For 
forced  vibration  problems  where  the  excitation  is  at  a  discrete  frequency,  the  solution  is 
relatively  simple  using  a  complex  modulus  if  only  the  forced  vibration  is  of  interest.  If  the 
spectrum  of  the  forcing  frequency  covers  a  spectrum  of  frequencies  then,  again,  a  different 
complex  modulus  should  be  used  at  each  freqency  involved.  McTavish  and  Hughes  [M2J  set  up 
a  mathematical  formulation  by  which  the  modulus  of  viscoelastic  material  is  represented  in  terms 
of  a  power  series  in  the  frequency  domain.  The  terms  of  the  series  are  factored  into  products  of 
quadratics,  each  quadratic  appearing  like  that  of  single-degree-of-freedom  system,  but  the  idea  is 
much  more  sophisticated  than  just  connecting  extra  oscillators  to  repre.sent  damping. 

From  the  theory  of  linear  viscoelasticity  |C.31,  the  constitutive  equation  which  could 
represent  the  one-dimensional  behavior  of  a  material  is 


a(t)  =  G(t)  e(o)  +  G(t-i) 


r. 


d£(T) 

dt 


e(i)  di . 


(4.25) 


Here  G(t)  is  the  stress  response  due  to  a  unit  step  strain  input,  which  is  the  definition  of  the 
material  relaxation  function. 

Taking  the  Laplace  Transform  of  Eq.  (4.25),  the  result  is  expressed  symbolically  as 


a  (s)  =  s  G(s)  e  (s)  (4.26) 

where  s  G(s)  is  the  material  dissipation  function. 

The  series  presented  in  detail  in  |M2)  is  as  follows,  where  the  fraction  represents  the 
behavior  of  over-damped  mini-oscillators,  with  the  proper  choice  of  parameters: 


s  G(s)  = 


(s  h-2i;^v 


k  “k  > 


where  to.  >  0  and  C  >  1 .  for  overdamping  in  the  mini-oscillators. 


(4.27) 


The  complex  modulus  G*(0))  =  j  to  G(j(o)  =  G  (to)  +  j  G'  (to) 
Tj  is  defined  by 


The  material  loss  factor 


G  (to) 

^  (to)  =  ^ -  .  (4.28) 

G  (to) 

In  1M2|,  the  generalized  viscoelastic  finite  element  is  represented  as  follows.  The 
global  linear  matrix  second-order  system  consists  of  the  following,  where  n,  m,  and  /  are 
integers: 

finite  elements,  subscripted  M’ . 
nil  material  moduli,  subscripted  ’j’  . 

Ilj  mini-oscillator  terms,  subscripted  'k' . 
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The  stiffness,  mass  and  damping  matrices  for  the  ith  element  are  K®,  M®  ,  and  Cf  req)ectively. 
The  element  stiffness  matrix  is  represented  as  a  linear  sum  of  contributions  from  the  m.  moduli 


m. 


K®  y'  GP  K®. 

.  u  U 


(4.29) 


where  G?  is  the  equilibrium  value  of  the  yth  modulus.  Further,  it  is  suggested  that  the  Kf 
matrices  may  be  diagonalized,  using  the  eigenvector  matrices  R.^ .  by 


R'  K  R  =A 

u  u  u  y 

where  A.^  is  a  diagonal  matrix.  Further,  the  relationships  exi.st  such  that 

r!*^  R  =  I,  so  K  =  R  A 
ij  ij  ij  y  y  y 

The  mini-oscillator  representation  for  the yth  material  of  the  /th  element  is 


sG..(s)= 
y  y 


y 


I 


k=l  ’  **  “ijk*  “ijk*  J 

4.2.3. 1  Finite  element  with  two  minimscillaturs  to  represent  viscoelastic  damping. 

The  finite  element  with  two  mini-oscillators  to  represent  viscoelastic  damping  will 
have  the  following  differential  equations 


(4.30) 


(4.31) 


(4.32) 


a.G^ 

-1 


a,G‘ 


“2 


q 

Z. 


>  + 


0 


a,G"2C 


0). 


a,G°2C 


2  j^e 


to. 


q 

Z| 

< 

1 

, 

G‘’(1  +a|  -i-a^ 

)K®  a,G"K®  tt^G^’K® 

q 

f(t) 

a,G‘’  K® 

a,G‘’  K®  0 

* 

*1 

►  =  < 

0 

a.,G‘’  K® 

0  a^G"  K® 

,  ^2, 

0 

[lat  the  bold-face. 

lower  case  letters,  here  q. 

L 

r 

*2 

and  fl 

(4.33) 

and  f(t),  represent  column 


matrices,  or  vectors. 

For  the  /th  bar  undergoing  only  axial  deformation,  the  stiffness  matrix  is 


A.  G"  r  1  1 

K.  =G‘’K®  =-i— !- 
•  '  •  /.  -11 
1  L 


(4.34) 


AE 


McTavish  and  Hughes  1 M21  suggest  diagonalizing  this  matrix.  Eigenvalues  are  ^  =0 , 2  — — . 
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The  matrix  of  normalized  eigenvectors  is 


«=-  '  ' 

VT  L-'  'J 

r’’’  K'  R  =  A  “I 

(4.35a) 

and 

(4.35b) 

/  Lo  oj 

If  the  diagonalized  matrix  is  used,  then  a  coordinate  change  is  involved  and  the 
physical  picture  of  an  actual  mini-oscillator  connected  in  the  system  bi^omes  more  obscure. 
However,  using  the  diagonalized  matrix  makes  the  computation  slightly  more  efficient  To  make 
the  change  in  coordinates,  let  z  =R  w  and  pre-multiply  the  second  and  third  set  of  equations  by 

R^.  Then  Eqs.  (4.33)  become 


a,G" 


GCl 


q 

'"i 

W., 


0  a,G^2C 


'  K" 


0) 


1 


K"' 


0), 


q 

w.. 


G"  1 1+(  ttj  +  a2  )|  K"' 

ajG^K*^  R 

a2G"K‘'‘R  ' 

q 

f(t) 

ajG^R*^  K'"' 

ajG”  A 

0 

'^1 

•  =  ' 

0 

» 

a^G^R”^  K'-' 

0 

a2G‘’  A 

.  '"2, 

0 

(4.36) 


The  matrix 


T  . 


2  -2 
0  0 


4.2.3.2  Representation  and  optimization  of  viscoelastic  damping  in  the  ten-bar  truss. 


Viscoelastic  material  might  be  applied  as  a  thin  layer  on  the  sides  of  a  bar  of  square 
cross-section.  The  present  analysis  assumes  that  there  is  an  outer,  rigid  constraining  layer, 
which  is  anchored  to  ground,  which  is  not  entirely  realistic.  A  more  reasonable  assumption 
would  be  to  represent  the  constraining  layer  as  a  rigid  mass,  free  to  move  in  the  direction  parallel 
to  the  axis  of  the  bar.  The  configuration  might  then  be  called  a  dynamic  absorber,  with  the 
viscoelastic  material  furnishing  the  stiffness  and  damping  while  the  constraining  layer  is  the 
absorber  mass,  which  is  the  idea  pursued  by  Carabetta  (Cl  ] 

A  reasonable  approximation  to  the  shape  of  the  real  and  imaginary  parts  of  a  practical 
damping  material  is  given  by  a  two-term  representation,  based  on  the  following  values  of  the 
parameters. 


a,  =0,5 

C  ,  =  15.0 

(1)  J  =  60,(XX) 

=  1.5 

c  2=  10.0 

=  10,000 
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If  the  viscoelastic  material  is  applied  to  the  side  of  a  truss  bar  and  it  is  deforming 
primarily  in  shear,  then  the  area  A.  would  be  the  glue  area  between  the  material  and  the  bar  and 

the  length  /.  would  be  the  thickness  of  the  viscoelastic  layer,  possibly  in  the  range  of  0.02S  to 

0.100  inches.  In  Table  4.8,  the  eigenvalues  are  listed  for  the  ten-bar  truss  with  two 
mini-oscillators  on  each  of  bars  3  and  4.  The  mini-oscillator  parameters  are  as  shown  above. 
Since  bars  3  and  4  are  of  different  cross-sectional  area,  the  A.//,  for  the  viscoelastic  layer  were 

taken  in  the  ratio 


0^/3 


=  0.8209 


The  modal  damping  ratios  and  natural  frequencies  are  listed  in  Table  4.8  for  three 


different  choices  for  a  parameter  "par^" 


G?A7/-  -  .  . 

defined  as  par,  =  ■  -  -  =  10^  10^  and  10^  . 
^  3  18.42 


The 


structural  modes  occur  in  complex  conjugate  pairs,  as  usual.  The  fictitious  modes  are 
overdamped  and  the  eigenvalues  are  real  and  negative.  It  is  difficult  to  decide  which  pairs  of 
overdamped  modes  to  take  together  for  the  purpose  of  computing  equivalent  damping  ratios  and 
natural  frequencies  for  the  fictitious  modes,  but  if  the  modes  are  not  strongly  coupled,  the 
eigenvectors  aid  in  identifying  the  pairs.  Thus  the  pairs  were  combined  using  the  following 
relationships 


2 

X  X  =0)  and  X,  +  A.  ,  =  2  C  co 

r  r+ 1  r  r  r+ 1  ’  r  r 

and  the  values  of  ^  ^  and  oj  ^  are  calculated,  with  the  results  as  given  in  Table  4.8.  The 

frequencies  of  the  mini-oscillators  themselves  were  1(),(K)()  and  6(),()()()  rad/s  so  the  results  of  the 
pairings  seem  reasonable. 


Table  4.8  Values  of  ^  ^  and  using  overdamped  mini-oscillators. 


par3= 


G^A^//^ 

18.42 


l(f 


Kf 


10 


5 


to 

r 

0) 

r 

(0 

r 

O.OOO  99 

131.17 

0.008  91 

131.69 

0.021  62 

136.85 

O.OOO  47 

173.42 

0.004  41 

173.75 

0.019  37 

178.10 

0.006  26 

270.76 

0.053  80 

281.73 

0.078  1 1 

366.37 

0.000  02 

326.05 

0.000  22 

326.07 

0.000  46 

325.72 

0.001  36 

462.56 

0.014  19 

467.50 

0.194  19 

642.52 

0.000  91 

600.52 

0.007  22 

604.55 

0.058  28 

556.67 

0.000  13 

665.05 

0.000  57 

665.69 

0.002  02 

664.77 

0.010  37 

809.80 

0.079  10 

940.63 

0.141  36 

1953.15 

60.77 

59,237. 

65.57 

54,905. 

68.90 

52,252. 

10.01 

59.962. 

10.06 

59,623. 

10.62 

56,438. 

15.11 

9,930. 

16.70 

8,979. 

36.39 

4,118. 

15.05 

9,968. 

15.53 

9,660. 

23.28 

6,440. 
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For  low  modulus  damping  material,  the  results  show  that  modes  3,  5  and  8  experience 

the  highest  damping  ratios,  in  the  first  column,  for  par^  =  10'^.  This  is  in  line  with  the  sensitivities 

calculated  in  Table  4.3  for  viscous  damping.  When  the  effective  material  modulus  is  increased 
by  a  factor  of  10,  the  C  ^  increase  by  almost  the  same  factor  for  the  structural  modes,  as  shown  in 

the  third  column  with  par.^  =  10^.  But  when  there  is  another  increase,  to  10^,  the  effect  becomes 
more  nonlinear  and  there  are  large  changes  in  the  to  ^’s.  There  is  a  stiffening  effect  from  the 
viscoelastic  material,  so  some  of  the  natural  frequencies  increase  while  others  decrease  as  the 
damping  layer  is  stiffened.  As  the  parameter  par.^  is  increased  to  10^,  the  modal  damping  ratios 
begin  to  decrease  so  those  shown  in  the  table  are  near  the  optimum  values. 


4.2.4  Design  and  uptimi/.atiun  of  viscoelastic  damping. 

The  use  of  added  viscoelastic  damping  is  desirable  but  the  behavior  is  still  hard  to 
predict,  for  many  reasons.  Testing  and  test  data  is  essential  and  it  is  difficult  to  perform  tests  of 
viscoelastic  material,  taking  into  account  the  preload,  temperature,  frequency,  previous  loading 
history,  and  damping  mechanism.  Having  the  material  properties,  then  the  following  factors 
must  be  considered.. 

1.  The  damping  mechanism.  The  vi.scoelastic  material  damps  best  when  deformed  in  shear. 

2.  Optimum  configuration.  In  a  truss,  load-carrying  members  could  include  a  slip  joint 
where  the  damping  material  transmits  the  axial  load  in  shear.  If  the  damping  material  has  to 
tran.smit  all  the  structural  load,  then  it  mu.st  have  appreciable  shear  strength  at  low 
frequencies.  Viscoelastic  damping  layers  might  be  glued  to  the  sides  of  bars  used  as  truss 
members,  with  a  constraining  layer  used  as  a  cover.  If  the  constraining  layer  has  the  same 
axial  strain  as  that  of  the  bar,  then  there  is  no  shear  deformation  of  the  viscoelastic  layer  and 
little  damping  is  produced.  Therefore  patches  of  viscoelastic  material  have  been  suggested 
by  some  investigators.  These  patches  then  have  the  appearance  of  dynamic  absorbers, 
where  the  viscoelastic  material  supports  the  absorber  mass,  which  is  that  of  the  constraining 
layer. 

2.  The  temperature  range.  Viscoelastic  materials  have  an  optimum  operating  temperature, 
so  material  must  be  chosen  which  has  a  high  loss  factor  in  that  range. 

3.  Heating  during  operation.  The  energy  is  dissipated  as  heat,  so  arrangements  must  be 
made  to  assure  that  heat  generated  can  vented  away. 

4.  Heating  during  fabrication.  Some  structural  fabrication  proce.sses  involve  heating.  Care 
must  be  taken  not  to  overheat  the  viscoela.stic  material.  For  example,  some  materials  are 

irreparably  damaged  if  heated  above  350*’  F. 


5.  Frequency  range  of  operation.  It  is  important  that  damping  material  is  chosen  which  has 
desirable  properties  in  the  frequency  range  where  damping  is  required. 

6.  Optimum  location  in  the  structure.  This  topic  is  addressed  in  Chapter  2.  The  optimum 
location  for  viscoelastic  damping  may,  or  may  not  be,  the  same  as  the  optimum  location  for 
viscous  damping. 

7.  Optimum  amount  of  viscoelastic  material.  This  has  to  do  with  the  area  and  thickness  of 
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individual  pads  of  material. 

With  regard  to  the  present  investigation,  in  addition  to  the  work  summarized  above,  the 
optimum  design  of  constrained  layer  damping  in  beams  was  a  topic  of  parallel  study.  Given  a 
beam  material,  such  as  aluminum  or  steel,  then  one  must  choose  the  damping  material  of 
optimum  thickness  and  optimum  Young’s  modulus.  Optimum  designs  are  discussed  in  [N8], 
based  on  formulas  given  in  |  N 1 1. 

A  senior  engineering  student.  Smith,  did  an  honors  thesis  [S2]  dealing  with 
optimization  of  the  thickness  and  elastic  modulus  of  an  added  damping  and  constraining  layer  in 
a  beam  in  bending.  He  made  a  finite  element  model  and  worked  in  the  frequency  domain,  with 
the  assumption  that  the  loss  factor  of  the  viscoelastic  material  was  rj  =1.0,  which  is  a  typical 
value  for  some  practical  materials.  He  found  effective  loss  factors  as  high  as  0.38  for  team 
tending  modes,  which  correspond  to  modal  damping  ratio  ^  =  0.19.  He  also  found  modes  in 
which  the  constraining  layer  moved  parallel  to  the  beam  axis,  which  was  very  efficient  for 
placing  the  damping  layer  in  shear,  with  maximum  ^  =  0.42.  An  unexpected  discovery  was 
that,  due  to  the  low  axial  stiffness  of  the  damping  layer,  additional  low  frequency  modes  are 
produced  in  which  the  constraining  layer  moves  out  of  phase  with  the  team.  The  maximum 
modal  damping  ratio  of  11-48  was  reported.  He  called  it  a  dynamic -absorber  type  mode, 
tecau.se  the  damping  layer  is  acting  like  an  absorber  spring,  but  it  is  a  stretching  of  the 
viscoelastic  material  which  allows  the  constraining  layer  to  move  perpendicular  to  the  axis  of  the 
team. 


In  another  phase  of  the  study,  a  senior  engineering  student,  Caratetta  [Cl],  studied 
damping  produced  by  four  dynamic  absorbers  attached  to  an  individual  clamped-free  truss  bar 
undergoing  axial  deformation.  The  dissipation  can  be  maximized  by  tuning  the  absorber 
frequency.  An  advantage  of  this  arrangement  is  that  the  viscoelastic  material  does  not  have  to 
carry  the  static  load  in  the  truss  member.  The  viscoelastic  material  provided  the  damping  and 
spring  stiffness,  in  shear,  and  a  rigid  mass  was  attached.  The  dynamic  modulus  of  the 
viscoelastic  material  was  represented  by  the  mini-oscillator  approach  |M2]  described  above, 
using  only  one  term  in  the  series.  Carabetta  predicted  modal  damping  ratios  of  C  ranging  from 
O.IO  to  0.57  Further  study  of  the  representation  of  practical  viscoelastic  material  properties  is 
needed,  but  the  results  of  this  study  of  the  dynamic  absorber  arrangement  are  very  encouraging. 
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5.  ACTIVE  DAMPING  AND  CONTROL. 


In  many  practical  situations,  it  is  advantageous  to  augment  passive  damping  with  a 
sufficient  level  of  active  control  to  achieve  specified  goals  for  the  eigenvalues.  If  the  purpose  of 
active  control  is  to  increase  the  damping  ratio  in  certain  modes,  then  the  poles  must  be  shifted. 
An  example  showing  the  effect,  on  eigenvalue  sensitivities,  of  superimposing  of  control  forces 
on  viscous  damping  forces  is  given  in  reference  1N9)  and  in  section  3.1.2. 

This  chapter  deals  with  the  design  of  passive  damping  and  active  control,  without  the 
use  of  observers  to  reconstruct  the  state.  Emphasis  is  on  designing  systems  with  desirable  modal 

damping,  using  feedback  control  of  the  form  -  G  i  -  H  z  .  The  control  may  devised  in  system 
space,  system  modal  space,  state  space  or  state  modal  space. 

Direct  output  feedback  control  in  system  space  is  explained  by  Meirovitch  [M4,  p.252].  It 
is  advantageous  to  have  as  many  actuators  as  sensors.  The  motion  is  sensed  and  control  forces 
proportional  to  velocity  and/or  displacement  are  fed  directly  back  to  the  structure.  For 
semidefinite  systems,  with  rigid  body  modes,  one  of  the  goals  may  be  to  stabilize  the  system 
with  displacement  feedback.  Quite  often,  the  reason  for  adding  control  forces  is  to  achieve  a 
desired  level  of  dissipation  for  each  mode,  which  means  to  place  the  closed  loop  poles  in 
specified  locations  in  the  complex  plane. 

An  advantage  of  using  modal  coordinates  is  that  the  differential  equations  of  motion  are 
decoupled.  If  they  are  recoupled  due  to  the  nature  of  the  control  forces,  the  control  is  called 
coupled  modal  control.  Another  change  in  coordinates  may  be  made  to  decouple  the  modal 
equations.  If  the  control  forces  are  chosen  so  the  modal  differential  equations  remain 
uncoupled,  the  control  is  called  natural  control  and  the  associated  damping  is  proportional.  If  the 
passive  damping  is  initially  small  enough  to  be  neglected,  natural  control  is  readily  designed  and 
an  example  is  given  in  IM3|  where  the  control  forces  are  chosen  in  system  modal  space  so  the 
closed  loop  natural  frequencies  are  the  same  as  those  of  the  open  loop  structure. 

Much  of  control  theory  was  developed  in  state  space,  but  :'.ot  by  engineers  concerned 
with  structural  dynamics  problems.  For  example,  Moore  IM61  discussed  the  possibility  of 
designing  the  entire  eigenstructure  and  presented  related  necessary  and  sufficient  conditions  for 
prescribing  both  eigenvalues  and  eigenvectors.  This  idea  was  expanded  by  Andry,  et  al,  [Al],  for 
three  types  of  feedback,  namely:  full  state,  output  and  constrained  output.  Neither  paper  dealt 
with  the  special  form  of  the  state  equations  for  structural  vibrations.  However,  an  important 
concept  is  that  the  mode  shapes,  as  well  as  the  eigenvalues,  may  be  designed.  Inman  [13]  gave  a 
related  example,  but  in  .system  modal  .space  and  not  state  modal  space.  Starek  and  Inman  [S5] 
presented  a  solution  in  state  space  to  the  so-called  inverse  problem,  where  mass,  damping  and 
stiffness  matrices  are  constructed  based  on  changing  all,  or  part  of,  the  system  eigenvalues  and 
eigenvectors,  with  rigid  body  modes  included.  They  extended  the  formulation  to  include  the 
case  where  only  a  portion  of  the  eigenvalues  and  eigenvectors  are  known.  To  assure  that  the 
constructed  system  mass,  stiffness,  and  damping  matrices  would  be  symmetric  when  solving  the 
inverse  problem,  Starek,  Inman,  and  Kress  |S6j  started  with  symmetric  matrices  and  applied 
constraints  to  the  chosen  eigenstructure.  Juang,  Lim  and  Junkins  |J1 1  considered  the  closed  loop 
control  problem  with  either  state  or  output  feedback  and  used  a  QR  decomposition  to  generate 
the  orthonormal  basis  that  spans  admissible  eigenvector  .space  corresponding  to  each  assigned 
eigenvalue.  Neubert  |N1 1,N12|  considered  the  problem  of  combined  damping  and  control  and 
full  state  feedback  and  presented  two  approaches,  one  based  on  using  eigenstructure  sensitivities 
and  the  other  based  on  a  perturbation  approach,  which  is  shown  to  be  related  to  the  inverse 
eigenvalue  problem. 

If  the  dimension  of  the  control  vector  is  le.ss  than  that  of  the  state  vector,  natural  control  is 
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not  possible  and  the  modes  are  coupled  by  the  control  forces.  If  only  a  limited  number  of  modes 
are  controlled,  then  the  question  of  spillover  ari.ses. 

In  modern  control  theory,  optimum  timewise  control  forces  may  be  found  by  way  of  an 
objective  function  and  solution  of  the  Riccati  equation.  In  the  present  chapter,  an  example  from 
[M4]  is  discussed  where  the  optimum  timewise  control  forces  are  determined  in  state  modal 
space.  These  control  forces  can  then  be  converted  to  system  space. 

The  contributions  of  the  present  study  that  are  summarized  in  this  chapter  are; 

1.  The  equations  for  natural  modal  control  in  system  space  are  given  and  related  to  the 

criteria  for  proportional  damping. 

2.  Examples  are  presented  for  designing  symmetric  or  unsymmetric  control  matrices 

using  the  eigenvalue  sensitivities. 

3.  One-step  perturbation  methods  are  presented  for  designing  the  control  and  damping 

matrices  of  a  structure  to  produce  desired  changes  in  the  eigenvalues  or  eigenvectors. 

4.  An  example  from  reference  (M4],  where  the  optimum  timewise  control  forces  were 

found  in  state  modal  space  by  way  of  the  Riccati  equation,  is  discussed  in  light 
of  the  present  work. 

5.1  Constant  gain,  closed  loop  control. 

5.1.1  Uncoupling  of  the  differential  equations  and  modal  control. 

In  system  space,  the  differential  equations  are  of  second  order.  Excluding  viscous 
damping,  Eqs.  (2. 1 ),  are  uncoupled  by  the  coordinate  transformation 

z(t)=<|)q(t)  (5.1) 

Substituting  Eq.  (5. 1 ),  pre-multiplying  Eq.  (2.1)  by  <I>  ,  and  using  <1>  1V1<1>  =  [MJ 

q  +  o'^'kO  q  =0”^  F(t)  =  N(t)  (5.2) 

lMnl^q  +  |MJ^  (w^l^q  =-[glq-lhlq  =  N(t)  (5.3) 

The  diagonal  matrices  on  the  left  side  of  Eq.  (5.3)  indicate  decoupling,  but  if  [g]  and  [h] 
are  non-diagonal,  the  equations  are  still  coupled.  If  |gl  and  (h)  are  diagonal,  which  is  called 
natural  control,  then  the  equations  are  uncoupled. 

If  the  control  matrices  in  N(t)  are  specified,  and  d>^  F(t)  =  N(t),  then  F(t)  is  found 
through  the  use  of  the  orthogonality  reiationship.s,  so 

F(t)  =  <!>■'’'  N(t)  =  MO  lM^r°  N(t)  =  -  MO  I M  [l  g]  q  +lhi  qj  (5.4) 
Changing  back  to  the  /.(t)  coordinates 

F(t)  =o  '^N(t)=  MOl  M  N(t)=-M  Ol  M  [l  gj  0'‘  i  +lhl  0‘'z]  (5.5) 

The  differential  equations  with  feedback  are  now 

Mz  +  MO|  M^|'^[|  g|0'U+lhI0‘‘7.!  +  Kz  =0. 

Mz  +  Gz  +|H  +  Kj  z  =  ( 


which  is 
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(5.6) 

(5.7) 


where 


(5.8) 


U  =  M  chi  M  J'^l  gj  <h''  and  H=  M  <h{  MJ'^[  h]  <h*‘ 

Pre-multiplying  Eq.  (5.6)  by  M  * 

>:+ch|  gl<h'‘i+lh)<h  ‘z)+M  ‘ Ki  =0.  (5.9) 

or  z  +  m'Gz+M*'|H  +  KU  =0.  (5.10) 

which  is  of  the  form 

Iz  +  C,  z  +  K,z  =0.  (5.11) 

S.U  Criterion  fur  proportional  damping  or  proportional  control. 

The  question  is  when  do  C I  and  K|  represent  proportional  damping  or  control. 

•} 

In  matrix  form  the  question  is  Cj  Kj  =  Kj  Cj  (5.12a) 

or  m  'gIVI''  (H  +  K1  =  M*' IH+K1M'‘g  (5.12b) 

Since  G  and  H  can  be  chosen  separately,  the  Eq.  (5.12b)  is  separated  into  two  parts: 

M''g  M*‘  H  I  M'*  H  M*‘  G  (5.13a) 

9 

M''gM‘‘k  =  M-'kM''g  (5.13b) 

Now  substitute  for  M''g  and  M'*  K  by  comparing  Eqs.  (5.10)  and  (5.9) 

<1>1  gich*‘ <h  |a)“|°<h''  =  <h  1  (O^l m  ' Mchl  MJ‘^{  g]<h*‘  (5.14a) 
which  reduces  to 

I  M  I  g|  [  0)^1  I  I  (0^1  ^  I  I  g|  (5.14b) 

The  answer  is  that  Eq.  (5.14b  _  not  represent  an  equality.  However,  the  two  sides  of 

the  equation  will  be  equal  if  |gj  is  du  '  ,  which  is  then  a  requirement  for  proportional  control. 

Then  Eq.  (5.14b)  becomes  the  following  equality 

I  M I  gl^  I  0)‘|  ^  =  I  o>;|  ^  I  M I  g|^  (5.15) 

The  advantage  of  proportional  damping  or  proportional  control  is  that  the  eigenvectors  <h  do  not 
change,  which  means  that  the  coordinate  transformations  which  decouple  the  original  matrices 
will  also  decouple  the  closed  loop  equations.  A  direct  way  to  check  that  this  is  true  is  to  simply 

T  T 

look  at  the  products  dJ  (i  <J>  and  OHO. 

The  latter  product  O^HO=  O^  M  0|  h|  O’*  O  (5.16) 

reduces  to  O^HO=|hl,  (5.17) 

and,  similarly,  0^(J0=|gl.  (5.18) 

So,  the  open  loop  eigenvectors  decouple  the  closed  loop  equations,  but  only  if  [h]  and  [g]  are 
diagonal  matrices. 
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5.1  J  One>step  methods  for  design  of  viscous  damping  and  control  in  state  space  by 
choosing  eigenvalues  and/or  eigenvectors. 

Natural  modal  control  and  the  closely  related  IMSC  method  1M4]  are  usually  presented 
as  a  one-step  method  for  designing  control  matrices  in  system  space  starting  with  plant  matrices 
where  there  is  no  viscous  damping  or  where  the  viscous  damping  is  proportion^.  Chie  of  the 
usual  features  is  that  there  is  no  change  in  mode  shape  from  those  the  undaii^)ed,  open  loop 
system.  A  perturbation  method  is  presented  in  IM5|  by  which  the  changes  are  predicted  in 
eigenvalues  and  eigenvectors  due  to  a  small  change  in  a  real  matrix  A.  The  changes  are  predicted 
using  a  eigenvectors  of  the  unperturbed  system.  A  different,  more  direct,  approach  was 
developed  as  part  of  the  present  study  and  are  presented  as  three  separate  methods. 

The  three  methods  require  only  one  step.  The  goal  is,  given  an  initial  system  with  or 
without  viscous  damping,  what  is  the  change  dA  in  the  matrix  A  to  achieve  the  desired 
eigenvalues  and  eigenvectors.  The  matrix  A  is  in  state  space  and  is  unsymmetric.  The  desired 
right  and  left  eigenvectors  U  and  V  and  the  diagonal  matrix  of  distinct  eigenvalues.  A,  are 
known  or  chosen. 

More  specifically,  given  the  desired  changes  dA  and  dU,  in  eigenvalues  and  eigenvectors, 
respectively,  find  the  required  change  dA  in  the  state  space  matrix  A.  The  three  methods  are 
arbitrarily  numbered  I,  II,  and  III.  They  are  as  follows: 

Method  I.  The  sum  |C  +  (j|,  initially  zero,  is  designed  to  achieve  a  desired  change  in 
eigenvalues  using  the  eigenvectors  and  mode  shapes  of  the  undamped  system. 

Method  II.  The  matrices  C I  and  of  the  closed  loop  system  are  designed  to  achieve 

desired  changes  in  the  eigenvalues  with  no  change  in  mode  shapes,  but  with  the  eigen¬ 
vectors  "corrected". 

Method  111.  The  matrices  Cj  and  Kj  of  the  closed  loop  system  are  designed  to  achieve 
desired  changes  in  the  eigenvalues,  mode  shapes  and  eigenvectors. 


First  of  all,  the  A  matrix  can  be  recovered  from  the  given  eigenproperties  from  the 
relationship 


UAV’''=A  (5.19) 

The  Eq.  (5.19)  is  not  new,  with  regard  to  recovering  the  original  A  from  the  eigen¬ 
properties  on  the  left  side.  However,  some  interest  has  been  apparent  lately  in  extending  the 
equation  to  specify  desired  eigenpropertie.s.  The  in  references  |S5J,  156]  and  [Jl]  are  especially 
relevant,  as  summarized  on  page  66.  In  the  pre.sent  work,  an  iteration  approach,  and  a  related 
FORTRAN  computer  program  were  written  early  in  the  study.  The  idea  was  to  specify  a 
"'“sired  increment  dA  in  A,  estimate  the  corresponding  increments  in  U  and  V,  and  use  (5.19) 
to  find  the  corresponding  dA  .  It  was  found  that  the  iteration  usually  converged  in  one  cycle  but 
the  essential  form  of  the  A  matrix  was  not  always  preserved.  Further  investigation  led  to  the 
formulation  of  the  problem  as  a  one  step  perturbation  approach. 

The  basis  of  the  perturbation  approach  is  to  take  small  variations  are  taken  in  the  matrices 
in  the  Eq.  (5.19  ),  as  follows. 

(U  +  dU)(A+dA)(V’^+dV)''’  =  A+dA  (5.20) 

Expanding  Eq.  (5.20)  gives 
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U  Av'^+A  dv'''+  dU  AV'’*-KiUA  dv'*' 

U  dA  v’^'+UdA  dv'*'+  dU  dAV^+dUdA  dV^  =  A  +  dA  (5.21) 

Taking  advantage  of  Eq.  (5.19)  and  dropping  terms  involving  products  of  two  or  three  small 
variations,  but  keeping  dU  A  dV^  in  some  cases,  the  remaining  terms  are 

U  dA  V^+U  A  dv’^+  dU  A  +  dU  A  dv'''=  dA  (5.22) 


If  the  changes  in  eigenvectors  are  small  enough,  the  increments  in  the  eigenvectors  may 
be  negligible.  Otherwise,  some  method  is  needed  to  estimate  dV  and  dU,  or  an  iteration  process 
must  be  used.  If  the  dA  is  chosen  such  that  there  is  no  change  in  eigenvectors,  then  the  result  is 
straightforward  to  carry  out,  and  is  simply 


UdAV=dA  fordU  =  dV=0 


(5.23) 


At  this  stage  it  is  important  to  divide  the  eigenvector  matrices  into  their  upper  and  lower 
portions,  as  in  Section  2.2.2  ,  that  is 


U  = 


u 

u 

u. 


and  recall  that,  in  general,  A  ,  so  taking  variations 


(5.24) 


U^+dU^  =  ( U,+d  )  (A+  dA)  (5.25) 

dU^  =  U,  dA+  dU,  A  +  dU,  dA  (5.26) 

So,  if  the  mode  shapes  do  not  change  and  d  =0 ,  d  U^=s  d  A  *  0.  So,  even  if  the 

mode  shapes  do  not  change,  the  eigenvectors  in  state  space  will  change  if  the  eigenvalues 
change,  b^'ause  the  upper  n  a  2n  portion  of  the  U  matrix  changes.  Further,  under  the  usual 
normalization,  the  fact  that  the  lower  right  vectors  do  not  change  may  go  unrecognized  unless  a 
special  normalization  procedure  is  used,  in  this  section  then,  the  rth  right  eigenvector  is 
normalized  using  the  n  x  n  mass  matrix  M  by 


I  X II  n  X  n  II X  I  t  X  i 

(5-  M  I 

Then  U  is  found  from  U  =  D,  A  (5.28) 

u  u  / 

and  V  from  V^=  U  *  .  (5.29) 


While  the  eigenvector  matrices  are  square,  and  are  assumed  to  be  invertible  because  the  vectors 
are  independent,  their  upper  and  lower  portions  are  rectangular,  of  size  n  x  2n  ,  and  are  not 
invertible. 


5.1.3.1  Method  I:  Design  of  the  sum  of  the  matrices,  [C+G] . 

In  this  method,  the  damping  and  control  matrices  are  initially  zero,  so  C=  G  =  H  =  0  and 
the  desired  change  dA  in  the  eigenvalue  matrix  is  given.  The  2n  x  2n  matrix  dA  consists  of 
four  n  X  n  submatrices  as  shown  here 
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r  dA 

uu 

1 

Q. 

0  ■ 

dA, 

L  /u 

1 - 

e 

0 

Taking  dU=  0  and  dV=  0,  there  results  from  Eqs.  (5.23)  and  (5.24) 


dA  =  U  dA  V"*" 


U  dA  V"*"  U  dA  yj  ' 
u  u  u  / 

U,  dA  1},  dA  yJ 
I  o'  '  J 


(5.30) 


(5.31) 


T 

Now  it  happens  dA  ^  0  but  it  is  in  the  wrong  position  in  the  dA  matrix,  so  in  Method  1, 
the  approach  is  to  add  it  to  dA^^  and  set  dA^^  =  0.  In  fact,  it  is  easily  shown  that,  if  C  ^  =  0 , 


U,  dA  v7  =  U  dA  (5.32) 

/  /  u  u 

which  follows  from  Eqs.  (2.39b)  and  (2.41a)  that  give 

U  =  U,A  and  AV‘^=  vj  whenC ,=  0  (5.33) 

Hence,  in  Method  I.  the  desired  change  in  is 

dC,  =  U,dA  v7+U  dA  (5.34a) 

with  dK,  =0  (5.34b) 


In  order  to  further  understand  the  effect  of  perturbation,  it  is  of  interest  to  perturb  the 
other  related  equations,  since  the  Eq.  (5.21)  came  from  perturbing  the  equation  UA  V  . 

The  new  eigenvectors,  U  +  dU  and  V+  dV,  would  be  normalized  so  there  results 

(V'‘'-K1V'‘')(U-K1U)  =  I  (5.35) 

T 

Cancelling  V  U  =  I  and  neglecting  the  products  of  increments,  Eq.  (5.35)  becomes, 

dv’^U  +  v'^dU  =0  (5.36) 

Taking  variations  in  the  characteristic  equations  involving  the  right  eigenvectors  yields 

(A+dA)(U-KlU)  =  (U-kHJ)(A+  d  A)  (5.37a) 

A  dC  +  dA  U  =U  dA  +  dU  A  (5.37b) 

Similarly,  for  the  left  eigenvectors, 

a”^  dV  +  dA’’'  V  =  V  dA  +  dV  a  (5.38a) 

dv’’'A  +  v’^'dA  =  dA  V’’’  +  A  dV’’’  (5.38b) 

Now,  in  the  iteration  process,  it  is  noticed  that,  sometimes,  the  resulting  dA^^  =dA^  ^ ,  so 
the  question  is,  under  what  conditions  does  this  occur. 

iix2ii  2iix2ii  2nxn  n\2ii  2nx2n  2nxn 

U  dA  v’’’  U,  dA  vT  ?  (5.39) 

Three  cases  are  of  interest: 
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Case  1 :  When  C|=<).  the  equality  holds.  This  is  seen  from  the  relationships  of 
Section  2.2.2 

Case  2:  The  equality  holds  for  proportional  damping,  where  C|K|  =  Cj. 

Case  3:  The  equality  does  not  hold  when  the  increment  dA  =  yA.  The  four  quarters  of  the 
dA  matrix  are 


U  dAv’^  =  -dC,  UdAVT=-dK. 

u  u  1  u  /  I 

(5.40a) 

L,  dA  v"''  =  d  I  U,  dA  vj  =  0 

So  for  this  special  case,  the  products  will  be 

(5.40b) 

U  dAv'^=-7C,  U,dAV'^=Yl 

(5.41a) 

U^dAVj  =-yK,  U^dAVj'  =0 

(5.41b) 

T 

A  difficulty  occurs  here  because  dA  =y  I,  which  means  that  the  redundant 
equation  in  state  space  is  not  maintained. 

The  above  becomes  more  clear  by  the  following  example. 

Example  5.1  One-step  design  of  the  [C+G]  matrix  using  Method  1. 

For  this  example,  the  mass  matrix  M,  K  and  the  state  space  matrices  M*  and  K*  are 

[2  ()1 

M  =  K  = 

[()  2j 


2 

0 

0 

0 

0 

0 

300 

-50 

0 

2 

0 

0 

* 

0 

0 

-50 

400 

K  = 

0 

0 

-3(K) 

50 

300 

-50 

0 

0 

0 

0 

50 

-400 

-50 

400 

0 

0 

The  characteristic  matrix  A,  the  eigenvalues  A  and  the  right  and  left  eigenvector  matrices  are 


0 

0 

-150 

25 

(0 +1 1.8171  i)  0 

0 

0 

0 

0 

25 

-200 

0 

(0-  11.8171  i) 

0 

0 

A  = 

A  = 

1 

0 

0 

0 

0 

0  (0-t- 1 4.5036  i) 

0 

0 

1 

0 

0 

0 

0  0 

(0-14.5036  i)  J 

3(X)  -50 
-50  4(M) 
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7.7199  i  -7.7199  i 

-3.9247  i  3.9247  i 

-0.0553  i 

0.0553  i  0.0187  i 

-0.0187  i 

3.1977  i  -3.1977  i 

9.4750  i  -9.4750  i 

V  = 

-0.0229  i 

0.0229  i  -0.0450  i 

0.0450  i 

0.6533 

0.6533 

-0.2706  -0.2706 

▼  ^ 

0.6533 

0.6533  -0.2706 

-0.2706 

0.2706 

0.2706 

0.6533  0.6533 

0.2706 

0.2706  0.6533 

0.6533 

Because  of  the  special  choice  of  the  mass  matrix  and  the  normalization  procedure,  it  happens  that 
=U^ ,  but  this  is  not  usually  true. 

For  later  reference,  the  sensitivity  matrices  are  also  recorded  here.  The  S  _  =  S,  and  S  .=  S-  . 

Z  \  4  3 


s 


1 


0.4268  -  ().(MKK)  i 
0.1768 -0.(KKK)i 
0.0000  +5.0433  i 
0.0000  +  2.0890  i 


0.1768  +  0.0000  1 
0.0732  +  0.0000  i 
0.0000  +2.0890  i 
O.(XXK)  +  0.8653  i 


O.tKXX)- 0.0361  i 
0.0000  -  0.01501 
0.4268  -  0.0000  i 
0.1768  -  0.00001 


0.0000  -  0.0150  i 
0.0000  -  0.0062  i 
0.1768  +  0.0000i 
0.0732  +  0.0000  i 


0.0732  +  O.(KMX)  i 
-0.1768 +0.(HKH)i 
0.0000+  1.0620  i 
O.OOtK)  -  2.5639  i 


-0.1768  +  0.00001 
0.4268  -  O.OOtK)  i 
O.(KKK)  -2.5639  i 
O.OOtK)  +6.1898  i 


0.0000  -  0.0050  i 
0.0000  +  0.01221 
t).0732  -  O.tXXK)  i  - 
-0.1768  +  0.00001 


0.0000  +0.0122  i 
0.0000  -  0.0294  i 
0.1768  -  0.00001 
0.4268  -  0.0000  i 


The  goal  now  is  to  determine  the  A+  dA  matrix  which  will  produce  the  following 
specified  changes  dA  in  the  eigenvalues: 


d  A 


(-0.30  +  0.01)  0  1) 

0  (-0.30  +  0.0  i)  0 

0  0  (-0.50  +0.0  i) 


0 

0 

0 


0  0  0  (-0.50 +<).()  i) 


Using  Eq.  (5.21 ),  assuming  temporarily  that  the  changes  in  eigenvectors  may  be  neglected. 


U  dA  =  dA  = 


-0.3293  0.0707 
().()7t)7  -0.4707 
0  0 

0  0 


0  0 

0  0 

-0.3293  0.0707 
().()707  -0.4707 


The  problem  is  that  the  A+dA  will  not  produce  a  correct  characteristic  equation,  since  the 
lower  right  corner  of  A+dA  should  be  zero.  However,  the  trace  of  A+dA  will  be  the  same  if  the 
lower  corner  is  added  to  the  upper  corner  and  ♦he  lower  corner  is  made  zero.  So  the  modified  dA 
matrix  will  be  designated  dA^  and  the  matrix  A^  =  A  +  dA^  . 
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’  -0.6586 

0.1414 

0 

o' 

’  -0.6586 

0.1414 

-150 

25* 

0.1414 

-0.9414 

0 

0 

A^  = 

0.1414 

-0.9414 

25 

-200 

0 

0 

0 

0 

G 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

The  eigenvalues  of  are 


-().3(K)()  +  i  11.8133  0  0  0 

0  -().3(KK)- i  11.8133  0  0 

0  0  -0.5000  +  114.4950  0 

0  0  0  -0.5000  -i  14.4950. 

The  real  parts  of  the  eigenvalues  are  those  desired  and  the  natural  frequencies,  cOj  =  1 1.8171  and 


0)^=  14.5036,  are  the  same  as  those  of  the  undamped,  unperturbed  matrix  A.  The  mode  shapes, 
in  .  are  unchanged  as  indicated  by  the  dli  listed  next,  but  the  upper  halves  of  the  right 


eigenvector  have  changed. 


dU=  Uq-U  = 


-0.1960  -0.(M)25  i 
-0.0812  -0.(K)10i 
0 

0 


-0.1960  +0.0025  i 
-0.0812+0.00101 
0 
0 


0.1353  +0.0023  i 

-0.3299  -0.0056  i 
0 

0 


0.1353  -0.0023  i 

-0.3299  +0.0056  i 
0 

0 


-0.1960  +  7.71741 
-0.0812  +  3.19671 
0.6533 
0.2706 


-0.1960-7.71741 

-0.0812-3.19671 

0.6533 

0.2706 


0.1353  -3.9223  i 
-0.3266  +9.4693  i 
-0.2706 

0.6533 


0.1353  -3.9223  i 
-0.3266-  9.4693  i 
-0.2706 
0.6533 


Again,  for  later  reference,  one  of  the  sensitivity  matrices  is  recorded. 


0.4268  +  10.0108 
0.1768  +  10.0045 
().()0(M)  +  i  5.0449 
0.0000  +  i  2.0897 


0.1768  +  10.0045 
0.0732 +  i0.(K)  19 
0.0000  +  i  2.0897 
O.OOOO  +  i  0.8656 


O.(KMK)-  i  0.0361 
0.0(M)0-i  0.0150 
0.4268-10.0108 
0.1768-10.0045 


0.0000  -  i  0.0150 
0.(K)00  -  i  0.0062 
0.1768-10.0045 
0.0732  -  i  0.0019 


Example  5.2  Same  example  as  5.1,  but  with  NxN  matrices  and  natural  control. 

If  the  mode  shapes  do  not  change,  then  the  same  problem  could  be  done  by  2x2  matrices. 


M= 


The  eigenvalue  matrices  for 


m  'k  are 


3(K)  -50' 
-50  4(M)_ 


^  39.6447 
wz 


1  *-[-<^-^533  -0.2706' 
2IO.3553J  ^-[-0.2706  0.6533 
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where  the  O  have  been  normalized  by 


M  <l>=  IVf^=I. 

If  the  modal  control  matrix  |gi  is  specified,  then  the  control  matrix  [Gj  in  the  original  coordinate 
system  is  calculated  with  the  results  as  follows 


-0.3  Ol  p_r-0.6586  0.1414' 
0  -O.SJ  0.1414  -0.9414_ 


If  this  dA  is  added  to  A.  the  resulting  eigenvalues  are  A-k1A,  as  follows  and  the  eigenvectors  are 
unchanged. 


A  H 


(-0.3+11.81331)  0  0  0 

0  (-0.3-11.81331)  0  0 

0  0  (-0.5+14.4950  1)  0 

0  0  0  (-0.5-14.950  i) 


5.13.2  Method  II:  Design  of  the  entire  dA  for  initially  damped  or  undamped  systems. 


Recall  that,  in  general.  A  ,  so  taking  variations, 

U^+dUy  =  (  U^+d  U/ )  (A+  dA)  (5.42) 

dU  =  U,  dA+  dU,  A  +  d  U,  dA  (5.43) 

Mill 

The  Eq.  (5.43)  gives  a  general  relationship  for  the  change  in  the  upper  right  eigenvectors 
and  shows  that  even  if  the  mode  shapes  do  not  change  and  d  =0  ,  d  11,=  d  A  0  if  the 

eigenvalues  change.  So,  even  if  the  mode  shapes  do  not  change,  the  eigenvectors  in  state  space 
will  change  if  the  eigenvalues  change,  because  the  upper  nx2n  portion  of  the  U  matrix  changes. 
Further,  under  the  usual  normalization,  the  fact  that  the  lower  right  vectors  do  not  change  may  go 
unrecognized  unless  a  special  normalization  procedure  is  used.  In  this  section,  the  right 
eigenvectors  are  normalized  using  the  n  x  n  mass  matrix  M  by,  for  the  \th  vector, 

/  .r  n  n x n  nx  /  I  x  t 

uj,  M  1=1  (5.44) 

Then  U  is  found  from  U  =  U,  A  (5.45) 

and  V  from  V^=-  U'*  .  (5.46) 

While  the  eigenvector  matrices  are  square,  and  are  assumed  to  be  invertible  because  the  vectors 
are  independent,  their  upper  and  lower  portions  are  rectangular,  of  size  nx2n,  and  are  not 
invertible. 

If  the  mode  shapes  do  not  change,  then  the  following  process  may  be  used,  which  does 
not  appear  in  the  literature.  Let  A^=A+dA  and  A^  =  A  +dA.  Then,  if  the  lower  right  vectors 

do  not  change. 
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U,[A-kIA] 

u, 


and 


v'^=U-‘ 

G  ^G 


(5.47) 


and  the  corresponding  state  matrix  for  the  closed  loop  system  is 

=  (5.48) 

where  dA  is  the  desired  shift  in  the  eigenvalues  and  dA  represents  the  required  change  in  A. 
Since  only  M'  'C^and  m  'Kq  are  needed,  it  is  more  efficient  to  use 


(5.49a) 
(5.49b) 

An  important  point  is  that  the  essential  form  of  the  A  matrix  is  maintained,  which  requires  that 

,T 


(5.50a) 
(5.50b) 

The  equations  (5.50)  are  easily  proved  by  multiplying  {2.41a)  and  (2.41b)  by  and  taking 
advantage  of  Eq.  (2.47). 


'''1'g,^g''g, 


Example  5.3  Design  of  damping  and  control  matrices  using  Method  U. 


For  this  example,  the  mass  matrix  M,  C,  and  K  matrices  are 


1  0 

M  = 

0  2 


C  = 


8 

-4 


K  = 


A  = 


300  -1(K) 

-100  400 

fie  right  and  left  eigenvector  matrices  of  A  are 


-8 

4  -300  100 

(-0.2656-1-12.8083  1)  0  0 

0 

2 

-1.5 

50  -150 

A  = 

0 

(-0.2656-  12.8083  1)  0 

0 

1 

0 

0  0 

0 

0  (-4..4844-t-17.7483i) 

0 

0 

1 

0  0 

0 

0  0  (-4..4844-17.7483i)  1 

-0.5288-5.5239 

i 

-4.6636-1-16.0712  i 

D  = 

0.3861-8.1868 

(conj)  0.4283-  5.4188  i  (conj) 

-0.4302-K).0502 

i 

0.9136-K).0319i 

-0.6395- 

0.0169  i 

-0.2927-1-  0.0498  i 

The  second  and  fourth  columns  of  U  are  the  complex  conjugates  of  the  first  and  third  columns, 
respectively,  and  are  denoted  by  the  symbol  "(conj)". 
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5.1.3.3  Method  III:  Design  uf  damping  and  cuntrui  matrices  tu  achieve  desired 
eigenvalues  and  eigenvectors. 

The  state  space  equations  whh  output  feedlyack  are 

I  net)  -  A„  n  (t)  =  B  u(t)  (5.51) 

where  B,  of  size  2n  x  p,  is  the  control  matrix  and  u(t)  is  the  p  \  I  vector  of  timewise  control 
inputs.  It  is  often  economical  to  have  the  dimension  of  u(t)  much  less  than  2n,  to  minimize  the 
number  of  controllers  needed.  The  output  w(t),  with  D  the  p  r  2n  output  matrix,  is 

w(t)  =  Dmt)  (5.52) 

If  u(t)  =  w(t).  Eq.  (5.51 )  becomes 
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I  n  (t)  -  Ao  T1  (t)  =  BD  net)  (5.53) 

The  Eqs.  (5.53)  and  (2.28)  have  a  similar  appearance,  but  there  may  be  more  freedom  in 
choosing  the  elements  of  the  two  matrices  B  and  0  than  the  single  matrix  R. 

Moore  [4J  is  given  credit  for  being  the  first  to  show  that  state  feedback  provides  freedom 
to  choose  not  only  eigenvalues,  but  eigenvectors.  He  derived  necessary  and  sufficient  conditions 
for  the  existence  of  a  D  which  yields  prescribed  eigenvalues  and  eigenvectors  for  the  closed  loop 
system  and  presented  a  procedure  for  computation  of  D.  The  theorem  states  that: 

Let  {X.  ,  /  =  /  -2n  )  be  a  self-conjugate  set  of  distinct  complex  numbers.  There  exists  a 
real  2n  x  m  matrix  D  such  that 

I  Ao+  BI)|  Vj  =  X.  V.  i  =  1, 2  . . .,  2n  (5.54) 

if,  and  only  if,  for  each  i. 

(i)  V.  are  a  linear  independent  set  of  vectors  in  the  space  of  complex  2n-vectors. 

(ii)  V.  =  V.  when  X.  =  X. . 

II  II 

(iii)  V.  €  span 

i 

If  D  exists  and  B  is  of  rank  p,  then  D  is  unique.  The  N;^  is  identified  from  the 

i 

definitions 


where  the  columns  of 


=  (XI  -  B  I  and 


partitioned  to  be  compatible  with 


for  the  product 


(5.55) 

,  furnish  a 


basis  for  the  null.space  of  . 


The  work  of  Srinathkumar  |8|  is  also  relative,  because  it  was  shown  in  [8]  that  for 
controllable  systems,  a  maximum  of  2nxp  eigenvector  entries  can  be  arbitrarily  specified  and  that 
no  more  than  p  entries  in  any  one  eigenvector  can  be  chosen  arbitrarily. 

In  the  present  paper,  we  choo.se  p  =  2n,  with  full  state  feed-back,  which  means  all  the 
eigenvalues  and  eigenvectors  may  be  specified,  but  care  must  be  taken  to  choose  independent 
eigenvectors  and  the  relationship  between  the  upper  and  lower  portions  of  the  eigenvectors  must 
be  kept  in  mind. 

In  general,  the  mode  shapes  and  eigenvalues  may  be  chosen,  with  the  restrictions  as 
outlined  above.  In  full  state  feedback,  all  the  entries  in  may  be  chosen  and  we  restrict  upper 

so  that  =11^  A.  The  operation  is  represented  by  Eqs.  (5.48)-(5.5()). 


Example  5.4  Design  of  damping  and  control  matrices  using  Method  Ill. 


For  this  example,  the  mass  matrix  M,  C,  and  K  matrices  are 


1  O' 
0  2 


C  = 


5  -4 
-4  3 


400  -200 
-100  700 


The  characteristic  matrix  A,  the  eigenvalues  A  and  the  right  and  left  eigenvector  matrices  of  A  are 


U  = 


-5 

4 

-4(K) 

2(X) 

(-1.4(K)4+18.6658ii)  0  0 

0 

2 

-1.5 

50 

-350 

0 

(-1.4004-18.66581)  0 

0 

A  = 

A  = 

1 

0 

0 

0 

0 

0  (-2.5996+27.116511) 

0 

0 
^  r 

1 

0 

0 

0 

0  0  (-2.5996-27.11651) 

-0.419+17.523  i  -0.4 1 9- 1 7  .523  i  5.798-10.0771  5.798+10.0771 
-1.937+  4.851  i  1.937-  4.851  i  -0.332+18.1771  -0.332-18.177  i 
0.935-  0.048  i  0.935  +  0.048  i  -0.389-  0.177  i  -0.389+  0.177  i 
0.266+  0.084  i  0.266-  0.084  i  0.665-  0.052  i  0.665+  0.052  i 


The  goal  is  to  design  the  mode  shapes  .so  the  normalized  modal  amplitudes  are 

.8165  0.8165  -0.5774  -0.5774] 

.4082  0.4082  0.5774  0.5774] 


u 

01  [()., 


in  order  to  make  the  amplitudes  in  the  ratio  of  2/1  for  the  first  two  modes  and  -1/1  for  the  last 
two.  The  desired  A+dA  is 


A+dA=^ 


(-6.0(K)0+18.6658i  i)  0  0  0 

0  (-6.(KXK)- 18.66581)  0  0 

0  0  (-8.00(M)+27.1165  i)  0 


0 


0 


0 


(-8.(KKK)-27. 11651) 


The  complete,  designed,  matrix  is  then 


and 


-4.899+ 1 5.24 1  i  -4.899- 1 5.24 1  i  4.6 1 9- 1 5.656  i  4.6 1 9+ 1 5.656  i 
-2.450+7.621  i  -2.450-  7.620  i  -4.619+15.6561  -4.619-15.656  i 
0.8165  0.8165  -0.5774  -0.5774 

0.4082  0.4082  0.5774  0.5774 

-13.33  2.67  -522.71  276.60 

1.33  -14.67  138.30  -661.01 
1  0  0  0 

0  1  0  0 


^G  '^G 


The  resulting  matrices  are 


-MC  =  C+dl  C 


-MK,=  K+d 


K  +  Hl=j^ 


r  13.33 

-  2.67] 

[  -  2.67 

14.67] 

522.7! 

276.60 

-276.60 

1322.02 
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5.1.4  Shifting  poles  using  eigenvalue  sensitivities,  symmetric  and  unsymmetric  control. 

Example  5.5  Design  of  symmetrical  control  matrix  using  sensitivities. 


Suppose  we  start  with  an  undamped  system  and  then  add  a  symmetrical  C+G  matrix. 


The  eigenvalues  are 
The  eigenvectors  are 


\)l  '  ^01“ 


A 


o 


0  0-54 
0  0  2  -2 
10  0  0 
0  1  0  0 


±0.5463i  and  2.5Sii7i 


U 

o 


0.2816  i  -0.2816  i  -2.2183  i  2.2183  i 
0.3310  i  -0.3310  i  0.9436  i  -0.9436  i 
0.5155  0.5155  -0.8569  -0.8569 

0.6059  0.6059  0.3645  0.3645 


The  sensitivities  and  are  given.  Note  that  8^2=  and  “  ^03 


0.1329 

0.1562  -0.2432  i 

-0.2859  i 

0.3123 

0.3671  0.5718  i 

-0.6720  i 

0.0726  i 

0.0853  i  0.1329 

0.1562 

0.1706  i 

0.2006  i  0.3123 

0.3671 

'  0.3671 

-0.1562  -0.1418  i 

0.0603  i ' 

-0.3123 

0.1329  0.1207  i 

-0.05131 

0.9504  i 

-  0.4043  i  0.3671 

-0.1562 

-0.8086  i 

0.34401-0.3123 

0.1329 

Now,  the  eigenvalues  desired  are 

±0.5399i  and  ^,,2  =-(>1680 ± 2.5842i 

The  changes  in  eigenvalues  are 

AXjjj  ,AX  -().()82()  +  ().{)()64  i  and  AX,j2  ,  ='(>•  l^l^O  +  0.0045  i 

The  complex  conjugates  provide  no  additional  information  with  regard  to  choice  of  changes  in 
the  plant  matrices.  So  there  are  four  knowns,  the  real  and  imaginary  parts  of  the  changes  in  two 
eigenvalues.  The  question  is,  which  four  matrix  elements  in  the  plant  matrices  should  be  chosen. 
We  choose  C(  1,1),  C(2,2).  K(l,l)and  K(2,2).  The  equations  are 

AX ,, ,  =  ( 1 , 1 )  AA(1 , 1 )  >  (2,2)  AA(2,2)  +  S^,,  ( 1 ,3)  AA(  1 ,3)  +Sq j{2,4)  AA(2,4) 

AX  =  S„3(l,l )  AA(1,1 )  +  Sj, 3(2,2)  AA(2,2)  +  S,^3(l,3)  AA(1,3)  +Sq3(2,4)  AA(2,4) 


0.1329  0.3671  -0.2432  i  -0.67201 

0.1329  0.3671  0.24321  0.67201 

0.3671  0.1329  -0.14181  -0.05131 
0.3671  0.1329  0.14181  0.05131 


T  f 

AA(1,1) 

,  AA(2,2)  . 
AA(1.3) 
AA(2.4) 


AA(1,1)  -0.4336 

AA{2,2)  _  -0.0664 

AA(1.3)  "  0.0326 

AA(2,4)  (-0.0023 

The  resulting  A  matrix  Is 

'-0.4336  0  -4.9674  4 

()  -0.0664  2  -2.0023 

I  1  0  0  0 

()  1  0  0  _ 

\)1  ’^01^  -0.0827  ±  0.53461  and  ±2.57341 

The  eigenvectors  are 

’-0.0367+0.27861  0.2317-2.20141 

^  -0.0527+0.32271  (conj)  0.041 8+0.9502.  (conj) 

'  0.5193-0.01171  -0.8577-0.03431 

0.6044+0.00501  0.3666-0.04011 

The  sensitivities  Sj  j  and  are  given.  Note  that  Sj2=S||  andSj^  =  Sj^ 

'0.1363+0.01481  0.1579+0.02211  -0.0114-0.25311  -0.0042-0.29481 ' 

g  ^  0.3159+0.04431  0.3658+0.06281  -0.0084-0.58961  0.0114-0.68611 

0.0399+0.07801  0.0436+0.09221  0.1313-0.09491  0.1561-0.10571 

-0.0288+0.16821  -0.0395+0.19461  0.3153+0.00511  0.3666+0.01731, 

'  0.3637+0.05311  -0.1579+0.0(8)61  0.0114-0.14211  0.(X)42+0.0611i  ' 

g  ^  -0.3159+0.00121  0.1342-0.02061  0.0084+0.12221  -0.0114-0.05141 

^13  -0.0399+0.95021  -0.0436-0.40631  0.3687-0.00851  -0.1561+0.02711 

0.0288-0.81311  0.0395+0.34731  -0.3153+0.(8)931  0.1334-0.02401, 
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Example  5.6  Design  uf  an  unsymmetricai  control  matrix  using  sensitivities. 

Given  the  Example  5.1  above,  suppose  it  is  desired  to  increase  the  damping  in 
eigenvalues  1  and  2  and  decrease  it  in  eigenvalues  3  and  4,  so  AX |=-  0.1  and  AX^  =  0.1  .  The 

starting  matrix  is  that  of  Example  5. 1  and  the  sensitivities  used  in  the  first  trial  arc  those  of 

the  undamped  system  The  only  element  to  be  changed  in  A  is  Aj2  •  From  the  sensitivity 

matrix  Sj  in  Example  5.1,  the  corresponding  sensitivity  is  0.1768.  To  find  AAj2  »  use  (0.1768) 

AA|2  =  -0.1  .T  he  result  is  AA^^  =  -  0.5656.  Make  this  change  in  the  matrix  A  and  call  the 

resulting  matrix  A^. 

-0.6586  -0.4242-150  25  * 

0.1414  -0.9414  25  -2(K) 

1  0  0  0 

0  1  0  0 

The  resulting  eigenvalues  are 

X,,  Xj  =  -0.3998  till  .8070  X  =-0.4(X)2  ±i  14.5022 

which  have  real  parts  very  close  to  those  desired.  If  an  iteration  process  is  used,  the  next  step 
would  be  to  update  the  sensitivities.  The  new  matrix  is 


S 


II 


0.4253  -(■  i  0.0262 
0. 1 787 -i  0.0293 
-0.2241  -t-i  5.0330 
0.6236-1-12.1048 


0.1755 -1-10.0177 
0.0744  -  i  0.(M)92 
-0.1741  -1-12.0779 
0.2238  -1-  i  0.8806 


0.0010-  i  0.0361 
-0.(M)30-i  0.0150 
0.4264  -I-  i  0.(K)45 
0.1763-10.0588 


0.0010 -i  0.0149 
-0.0010 -i  0.0063 
0.1763  +  10.0088 
0.0739-10.0215 


which  shows  that  the  sensitivities  are  now  complex  numbers,  so  a  corresponding  change  in  the 
imaginary  part  snould  be  demanded.  The  damping  is  not  proportionai  and  a  change  in  mode 
shape  is  expected.  The  new  right  eigenvectors  are: 


U 


I 


-0.2171  +17.7170  -0.2171  -17.7170 
-0.2147  +  i  3.12881  -0.2147  -  i  3.1881 
0.6535  -  i  ().()()37  0.6535  +  i  0.{K)37 
0.2703  +i  0.0090  0.2703  -  i  0.(M)90 


-0.8304  -  i  3.9956  -0.8304  +  i  3.9956 
-0.6548  +  i  9.4990  -0.6548  -  i  9.4990 
-0.2737  +  i  0.0648  -0.2737  -  i  0.0648 
0.6557  +  i  0.027 1  0.6557  -  i  0.027 1 


and  dlJ  =  U,-U„ 

I  o 


-0.0212 -i0.(H)04 
-0.1335 -i  0.0086  (conj) 
0.(HK)2  -  i  0.0037 
-0.(HM)3  +  i  0.0090 


-0.9567  -  i  0.0733 
-0.3282  +  i  0.0296  (conj) 
-0.(K)31  +10.0648 
0.(M)25  +  10.0271 
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Example  5.7  Same  problem  as  Example  5.6,  but  using  a  variation  of  Method  L 

For  this  example,  which  is  the  same  problem  as  done  in  Example  5.6,  the  idea  of  Method 
I  is  used,  but  the  eigenvectors  used  are  Uq  and  of  the  damped  system  of  Example  5.1.  The 

desired  eigenvalues  are  those  given  in  ,  namely 

11.8133  0  0  0 

-0.4- i  11.8133  0  0 

0  -0.4  +  114.4950  0 

0  0  -0.4  -  i  14.4950 


Aq  -I 


-0.4  +  i 

0 

0 

0 


-0.7293  0.0707 
0.0707  -0.8707 
1.0000  0.0000 
0.(MM)1  l.(KMX) 


-150.00  25.00 

25.00  -200.00 
-0.0707  -0.0707 
-0.0707  0.0707 


The  2  .r  2  matrix  in  the  lower  right  corner  of  AGj  is  non-zeio  and  not  equal  to  the  2  x  2  matrix 

in  the  upper  left  corner.  Thus,  in  the  spirit  of  Method  1,  we  add  it  to  the  upper  left  comer  and 
zero  the  lower  right  corner,  producing  the  modified  characteristic  matrix. 


Modified  AGj  = 


■  -0.8(MM)  0.00(K) 
O.(KKK)  -0.80(K) 
l.tKXX)  O.O(KM) 
O.O(KM)  l.(KHK) 


-150.00  25.00 

25.(M)  -200.00 
O.tKXK)  O.OOIX) 
O.CMKX)  0.0000 


5.U  Nonsymmetric  cuntrul  matrices,  forced  vibration  problem. 

To  do  a  forced  vibration  problem  with  symmetric,  or  unsymmeiric,  control  matrices  G 
and  H,  the  closed-loop  matrix  equations  are 


’  M  o' 

* 

i  (t)  1 

-  + 

'  IC+G! 

1K+H|1  fz(t)1 

<  »  =< 

f(t) 

(5.56) 

(1  -  K 

. 

K 

0  J  lz(t)J 

0 

The  shorthand  form  of  Eq.  (5.56)  is 

M*  j]  +K*  x\  =  F(t)  (5.57) 

or  n-A  n=|M*l''  F(t)  (5.58) 

where  A  is  the  characteristic  matrix  for  the  closed-loop  system  and  U,  V,  and  A  are. 
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respectively,  the  right  eigenvector,  left  eigenvector  and  eigenvalue  ntatrices  of  the  closed-loop 
system. 

The  equations  are  decoupled  by  substituting 

q  (t)  =  U  q(t)  (5.59) 

T 

into  Eq.  (5.58)  and  pre-multiplying  by  V  . 


5.2  Optimum  damping  and  cuntrul  by  way  of  performance  index  and  the  Riccati  equation. 

5.2.1  Performance  index  and  optimal  control 

An  often-used  form  of  a  quadratic  performance  index  for  a  controllable  system  isis 

J  =  i  I  z\t)  Q(t)  z(t)  +  z(t)  R(t)  z(t)l  dt  (5.61) 

where  Q  is  real,  symmetric  and  positive  semidefinite  while  R  is  real,  symmetric  and  positive 
definite.  Finding  the  optimal  solution  involves  solving  the  nonlinear  Riccati  equation,  but  which 
can  be  transformed  to  a  linear  form.  The  behavior  is  a  function  of  time,  but  at  t=t^  the  Riccati 

matrix  approaches  a  steady-state  value.  An  numerical  example  of  an  optimal  timewise  solution 
is  presented  on  page  199  in  reference  [M41. 

5.2.2  The  Riccati  equation  and  modal  control. 

If  the  system  is  undamped  and  the  control  is  such  that  the  xth  modal  equation  is  inde¬ 
pendent  of  the  other  modal  equations,  the  performance  index  may  be  written  for  each  mode  and 
the  total  performance  index  is  obtained  by  summing  the  modal  performance  indices.  For  the  Uh 
mode,  the  performance  index  is  given  by 


'r '  Jj[  “J*"  Qr  “r ^  “r  *'>] 
and  the  total  performance  index  J  is  the  sum  of  the  modal  performance  indices. 


(5.62) 


n 

J = X  ^ 

r=l 

The  can  be  chosen  so  the  first  term  in  Eq.  (1)  is  the  .sum  of  the  modal  kinetic  and  potential 
energies  of  the  plant.  Again,  the  state  space  differential  equations  are 

M*  n  +K*  T|  =  F(t)  =  G*  li  +  H*  T)  (5.64) 

Assume  the  eigenvalues  A  and  the  right  and  left  eigenvectors  LI  and  V  for  the  M  *  K  matrix  are 
available.  Hence  we  substitute  (5.65)  and  pre-multiply  (5.64)  by  M 

q  (t)  =  U  q(t)  (5.65) 

with  the  result 


U  q(t)  U  q(t)  =  M*'*  F(t) 


(5.66) 
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T 

To  decouple  the  equations,  pre-multiply  by  V  and  take  advantage  of  the  usual 


orthogonality  relationships. 

q(t)  -  Aq(t)=v'’'lVI*'‘F(t)=v'’'M*'‘GUq  +v'^M*’'HUq  (5.67) 

q  (t)  -  A  q(t)  =  Q(t)  =  M*‘‘  F(t)  =  |g|  q(t)  +lh|  q(t)  (5.68) 

In  general,  the  |g|  and  |h|  matrices  will  not  be  diagonal.  If  they  are  diagonal,  the  modal 
control  is  natural.  The  real  and  imaginary  parts  of  the  modal  control  force  Q(t)  are  designated 
X(t)  and  Y(t),  written  as 

Q(t)=X(t)  +  i  Y(t)  =  M*-'  F(t)  =  |g|  q(t)  +[hl  q(t)  (5.69) 

To  return  to  the  system  coordinates,  substitute  from  (5.65)  for  q(t) 

v'*'  M*'‘  F(t)  =(g|  U*'  i]  +|h)  U'*  r\  (5.70) 

which  gives  F(t)  =M  V'^^lgJ  U  *  f} +(h)  U  *  q  J  (5.71) 

or  alternatively  F(t)  =M  V‘^^X(t)  +  i  Y(t)J  (5.72) 


The  Example  5.X  is  done  next,  to  show  how  the  optimum  timewise  modal  control  forces 
and  associated  closed  loop  poles,  found  from  a  solution  of  the  Riccati  equation  are  converted  to 
control  forces  in  the  system  coordinate  system.  In  Example  5.9,  the  results  are  checked  by  the 
natural  modal  control  approach,  given  the  desired  clo.sed  loop  eigenvalues.  In  example  5.10,  the 
same  results  are  gotten  by  using  Method  11. 


Example  5.8  Conversion  of  modal  control  forces  to  system  control  forces. 

The  purpose  of  this  example  is  to  show  how  the  optimum  solution  obtained  from  the 
Riccati  equation  for  a  specific  problem  by  Meirovitch  |M4|  related  to  the  present  discussion. 
The  plant  matrices  are 

'‘“['4 1] 

and  the  eigenvalues  are  X|,  ±  i  0.5463  and  Xy  X.^=  ±  i  2.5887.  In  [M4]  special 

coordinates  are  introduced  so  that  all  the  matrices  will  contain  real  numbers,  rather  than 
imaginary  numbers.  The  (j  and  R  matrices  used  in  the  special  form  of  the  performance  index 
are  .aken  as  the  unit  matrir.  The  real  part  of  the  control  forces  X(t)  =0  and  the  imagl.iary  part  is 
of  the  form 

f^,(t) 

^11  ^12  ^  ]  I  ^1^*^ 

I 

Having  the  control  forces,  the  closed  loop  eigenvalues  can  be  determined.  An  alternate  approach 
to  that  used  in  |M4|  is  to  convert  the  timewise  control  forces  in  modal  coordinates  to  system 


coordinates.  Then  the  G  and  H  control  matrices  can  be  identified  directly.  Tlie  4,(0  and  (t) 
are  related  to  q(t)  by 


q,(t)  =  4,(0  +  «  n,(0  and  =  4,(0  -  i  n,(t) 

To  convert  back  to  complex  system  coordinates,  the  complex  conjugates  come  into  play.  The 
conversions  are  written  in  matrix  form,  like 


with  the  inverse 


0.5  0.5  1 

-0.5  i  -0.5  i  X 


The  control  forces  in  system  coordinates  are 


Q(t)  = 


-0.4404  0.0990 
0.0991  -0.7323 


z,(t) 

-1.3705  0.0738  1  z^O) 
0.0738  -2.6304  J  z^(t) 

Z2(t). 


and  the  G  and  H  control  matrices  are  identified  as 


G= 


0.4404  -  0.0990  ' 
-0.0991  0.7323 


and  H  = 


’  1.3705 
,-0.0738 


-  0.0738* 
2.6304 


From  this,  the  new  A  matrix  is  formed  and  the  closed  loop  eigenvalues  calculated. 
The  closed  loop  eigenvalues  are 


Xj  ,  X.2  = -0.6419  ±  i  0.4588  and 

so  <i)j=  0.5463  and 

The  changes  in  the  eigenvalues  are 

dX,= -0.6419  -  i  0.0875  and 
dX^  = -0.7009  -i  0.0017  and 


X-  ,  X  =  -0.7009  ±  i  2.5870 
3  4 

=  2.5887 

dX2  = -0.6419  -Hi  0.0875 
dX  =-0.7(K)9-h  i  0.0017 

4 
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Example  5.9  Same  problem  as  5.8,  but  using  pole  shifting  by  natural  control  method. 

The  goal  of  this  example  is  to  determine  the  required  control  matrices  using  the  method 
of  natural  modal  control  with  //  x  n  matrices,  rather  than  in  state  space.  The  form  of  the 
polynomial  in  A.  is 


gX+(h+(o")=0 
r  r  r  r 

and  the  given  roots  are 


or 

so 

and 


.X 


X  .X  =  a  ±  i 

r  r  r 


±  i 


(h  +Q} 
r  r 


-2  a 

r 


The  individual  eigenvectors  will  be  designated  as  and  the  //  x  n  eigenvector  matrix  as 
<I>. 


0.5155  0.8569 
0.6059  -0.3645 


The  required  diagonal  modal  control  matrices  are 

[1.2838  0  1 

^  [  0  1.4018J 

hJ  0.3241  0  1 

L  0  0.4825  J 


The  corresponding  control  matrices  in  the  original  coordinates  are 


1.3704  -0.0737 
-0.0737  2.6304 


and  G= 


0.4404  -0.0990 
-0.0990  0.7324 


so  the  K-t-H  matrix  is 


5.4404  -4.0990 
-2.0495  2.3662 
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Example  5.10  Same  as  but  using  Method  11. 

Suppose  we  start  with  an  undamped  system  and  then  add  a  symmetrical  C+G  matrix. 


0 

0 

-5 

4 

0 

0 

2 

-2 

1 

0 

0 

0 

0 

1 

0 

0 

li 

and 

^02’ 

The  eigenvalues  are  ±(l.5463i  and  ±  2.5887i 

The  eigenvectors  are 

’  0.2816  i  -0.2816  i  -2.2183  i  2.2183  i' 

0.3310  i  -0.3310  i  0.9436  i  -0.9436  i 

0.5155  0.5155  -0.85^9  -0.8569 

().6059  0.6059  0.3645  0.3645  . 

The  desired  change  in  A  is 


-0.6419-0.0875  i 


-0.6419  -K).0875 
0  -0. 
0 


0  0 

'5  i  0  0 

0.7009 -iO.tM)  17  0 

0  -0.7009 -i-i  0.0017 


-0.3309  0.2365  i  -0.6006  2.2168i 

-0.3889 0.2780  i  (conj.)  0.2555-0.094301  (conj.) 
0.5155  -0.8569 

0.6059  0.3645 


Ul  (A-KlAj 


^G  = 


-0.5618  i 
-1.3207  i 
0.2577  -  0.3606  i 
0.6059  -  0.8477  i 


-0.1656  i 

(conj.)  0. 1409i  (conj.) 
0.4284-  0.1161  i 
-0.3645  0.0988  i 


A^=U^|A-KdA|V‘ 


'^G  = 


-1.3704  0.0737  -5.4404  4.0989 

0.0369  -1.3152  2.0495  -2.3662 
1  0  0  0 

0  1  0  0 
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5.3  Summary  and  conclusions. 

The  present  study  has  covered  a  range  of  problems  related  to  the  efficient  representation 
and  design  of  damping  and  control  forces  for  discrete  systems.  Some  conclusions  are  given  at 
appropriate  locations  in  the  preceding  portions  of  the  report.  A  summary  is  given  here  in  the 
form  of  suggestions  for  further  study,  as  follows: 

1.  Some  forms  of  passive  and  active  damping  are  achieved  by  devices  representable  by 
viscous  forces.  However,  viscoelastic  damping  is  widely  used  and  more  work  is  needed 
to  develop  more  efficient  analytical  and  computational  techniques.  The  mini-oscillator 
approach  is  readily  adapted  to  present  modal  analysis  methods  and  computer  programs, 
but  a  considerable  increase  in  size  of  matrices  is  required.  The  fractional  calculus 
approach  may  be  more  efficient,  but  is  not  as  easy  to  use.  Complex  numbers  are  readily 
handled  in  present  computer  programs  and  the  use  of  a  complex  Young’s  modulus  is 
efficient  in  some  instances,  especially  if  the  modulus  is  fairly  constant  over  the  frequency 
range  involved. 

2.  Eigenvalue  and  eigenvector  sensitivities  may  be  used  for  constrained  optimization 
problems.  Their  use  to  achieve  certain  eigenvalues  through  design  of  dissipative 
damping  and  control  matrices  has  been  demonstrated  in  the  present  study.  They  can  also 
be  used  to  design  eigenvectors,  but  eigenvectors  tend  to  change  rapidly  with  changes  in 
parameters  and  more  study  is  needed  in  this  area. 

3.  Direct  design  of  damping  and  control  matrices  through  eigenvector  as  well  as  eigen¬ 
value  assignment  is  very  promising,  as  indicate  by  the  methods  and  examples  outlin^  in 
Chapter  5.  The  literature  on  the  inverse  eigenstructure  problem  is  relative  and  is 
increasing  rapidly.  Eigenvector  sensitivities  can  be  useful  in  choosing  the  eigenvector 
matrices  and  is  an  area  worth  further  development. 

4.  Practical  realization  of  the  system  and  control  matrices,  while  not  a  major  topic  of  the 
present  report,  has  been  mentioned  in  some  of  the  papers  published  by  the  writer  and  in 
many  other  papers.  Maintaining  data  bases  of  properties  of  viscoelastic  materials  in  forms 
that  are  readily  useful  in  computer  programs  for  structural  analysis  is  an  ongoing  task. 
Graphical  representation  is  helpful,  but  a  listing  of  the  parameters  involved  in 
curve-fitting,  as  described  in  Chapter  3,  would  be  directly  useful  to  the  analyst. 


% 


-90- 


REFERENCES 


Ai.  Andry,  A.  N.,  Jr.,  E.  Y.  Shapiro  and  J.  C.  Chung,  "Eigenstructure  Assignment  for  Linear 
Systems,"  IEEE  Transactions  on  Aerospace  and  Electronic  Systems,  AES- 19,  no.  5, 
pp.  713-729,  1983. 

Bl.  Bagley,  D.  L.  and  Torvik,  P.  J.,  "Fractional  Calculus-A  Different  Approach  to  the  Analysis 
of  Vicoelastically  Damped  Structures,”  AJAAJnl..  v.  21,  n.  5,  pp.  741-748,  1983. 

Cl .  Carabetta,  Christopher  J.,  "Damping  Analysis  of  Axial  Elements  Treated  with  Constrained, 
Sectional,  Viscoelastic  Layers,"  Senior  Honors  Thesis,  The  Penn.sylvania  State  University, 
May,  1991 

C2.  Caughey,  T.  K.  and  O'Keliy,  M.M.J.,  "Classical  Normal  Modes  of  Damped  Linear  Dynamic 
Systems,"  ./.  ofAppl.  Mech,  AS  ME,  Vol.  32,  pp.  583-588,  1965. 

C3.  Christensen,  R.  M.,  Theory  ofVisoelasticity:  An  Introduction,  2nd  Edition,  Academic  Press, 
New  York,  1982. 

FI.  von  Rotow,  A.  H.  and  D.W.  Vos,  "The  Need  for  Passive  Damping  in  Feedback  Controlled 
Rexible  Structures,"  Proceedings  of  Damping  '91,  ^Flight  Dynamics  Laboratory  of  the 
Air  Force  Wright  Aeronautical  Laboratories  1991)  pp.  GBB-1  to  GBB-12. 

Gl.  Gilheany,  J.  J.,  "Optimal  Selection  of  Dampers  for  Freely  Vibrating  Multidegree  of  Freedom 
Systems,"  In  Proceedings  of  Damping  '  89  (Flight  Dynamics  Laboratory  of  the  Air  Force 
Wright  Aeronautical  Laboratories  1989)  pp.  FCC-1  to  FCC-18. 

11.  Inman,  D.  J.  and  A,  N.  Andry,  Jr.,  "A  Procedure  for  Designing  Overdamped  Lumped 
Parameter  Sy.stem.s,"  Shock  and  Vibration  Bulletin,  54(5),  49-53  (1980). 

12.  Inman,  D.  J.  and  Andreyi,  A.  N.  Jr.,  "Some  Results  on  the  Nature  of  Eigenvalues  of  Dis¬ 
crete  Damped  Linear  Systems,"  Jnl.  of  Applied  Mechanics,  ASME,  47,  no.4,  pp.  927- 
930,  1980. 

13.  Inman,  D.,  Vibration  Control,  Measurement  and  Stability,  Prentice-Hall,  Englewood  Cliffs, 
NJ,  1989. 

Jl.  Juang,  J.,  Lim,  K.  B.,  and  Junkins,  John  L.,  1989,  "Robust  Eigensystem  Assignment  for 
Flexible  Structures,"  AIAA  Journal  of  Guidance,  Control  and  Dynamics,  Vol.  12,  No.  3, 
pp.  381-387. 

LI.  Liang,  Z.,  Tong,  M.,  and  Lee,  G.  C.,  "A  Strong  Criterion  for  Testing  Proportionally 

Damped  Systems,"  Proceedings  of  Damping  '91,  San  Diego,  Wright- Patterson  Air  Force 
Base  WL-TR-9 1-3078.  paper  ECD,  August  1991. 

M 1 .  McLoughlin.  F.  A.,  Integrated  Structural  Damping  and  Control  System  Design  for  High- 
Order  Flexible  Systems,  Ph.D.  Dissertation  (Stanford  University,  1991) 

M2.  McTavish,  D.  J.  and  Hughes,  P.  C.,  "Finite  Element  Modeling  of  Linear  Viscoelastic 
Structures,"  in  Rogers,  L.  and  Simonis,  J.  C.,  Editors,  The  Role  of  Damping  in  Vibration 
and  Noise  Control,  ASME  DE-Vol.  4,  Sept.,  1987. 

M3.  Meirovitch,  L.,  Introduction  to  Dynamics  and  Control,  iohnVi/'\\ey  &  Sons, 

New  York,  1985. 


-91- 


M4.  Meirovitch,  L..  Dynamics  and  Control  of  Structures,  Wiley,  New  York,  1990. 


M5.  Meirovitch,  L.,  Computational  Methods  in  Structural  Dynamics,  Sijthoff  &  Nordhoff, 

Rockville,  Maryland,  1980. 

M6.  Moore,  B.  C..  "On  the  Flexibility  Offered  by  State  Feedback  in  Multivariable  Systems 
Beyond  Closed  Loop  Eigenvalue  Assignment,"  IEEE  Transactions  on  Automatic  Control, 

AC-21,  pp.  689-692,  1976. 

Nl.  Na.shif,  A.D.,  et  al.  Vibration  Dampint;.  John  Wiley  &  Sons,  New  York,  1985. 

N2.  Nelson,  R.  B.,  "Simplified  Calculation  of  Eigenvector  Derivatives,"  AIAA  Jnl.,  v.  14,  n.9, 
pp.  1201-1205,  Sept.  1976. 

N3.  Neubert,  V.  H.,  Mechanical  Impdance.  ModelHn^l Analysis  of  Structures,  Jostens  Printing 
Company,  State  College,  PA,  1987. 

N4.  Neubert,  Vernon  H.,  "Optimization  of  Energy  Dissipation  Rate  in  Structures,"  Proceedings 
DAMPING  H9,  Flight  Dynamics  Laboratory  of  the  Air  Force  Wright  Aeronautical  Labora¬ 
tories,  1989,  pp.  FCD-I  to  FCD-24. 

N5.  Neubert,  Vernon  H.,  "Optimization  of  Location  and  Amount  of  Viscous  Damping  to  Mini¬ 
mize  Random  Vibrations,"  Oral  presentation  at  Acoustical  Society  of  America  Semi- 
Annual  Meeting,  May  1990. 

N6.  Neubert,  Vernon  H.,  "Optimization  of  Location  and  Amount  of  Viscous  Damping  to 
Minimize  Random  Vibrations,"  JASA,  v.  93,  n.  5,  pp.  2707-2715,  May  1993. 

N7.  Neubert,  Vernon  H.,  "Fitting  Curves  to  Data  for  Complex  Moduli  of  Viscoelastic  Materials, 

Technical  Report,  Applied  Research  Laboratory,  The  Pennsylvania  State  University,  (In 
preparation  from  a  completed  internal  memorandum.) 

N8.  Neubert,  Vernon  H.,  "Optimum  Design  of  Constrained  Layer  Viscoelastic  Damping  for 
Structures,"  Technical  Report,  Applied  Research  Laboratory,  The  Pennsylvania  State  Uni- 
sity,  (In  preparation  from  a  completed  internal  memorandum.) 

N9.  Neubert,  Vernon  H.,  "Sensitivity  of  Modal  Parameters  to  Nonsymmetrical,  Dissipative 
Control  Matrices,"  Proceedings  of  the  63rd  Shock  and  Vibration  Symposium,  October, 

1992,  Las  Cruces,  New  Mexico,  v.  1.  pp.  161-171.  >1 

NIO.  Neubert,  Vernon  H.,  "Sensitivity  of  Eigenvalues  to  Nonsymmetrical,  Dissipative  Control 

Matrices,"  Journal  of  Shock  and  Vibration,  v.  1,  n.  2,  Dec.,  1993  » 

Nil.  Neubert,  Vernon  H.,  "Design  of  Control  Matrices  in  the  Pre.sence  of  Viscous  Damping," 

Proceedings  of  the  64th  Shock  and  Vibration  Symposium.  Oct.,  1993,  Ft.  Walton  Beach, 

Ra.,  V.  1,  pp.  515-524. 


N12.  Neubert,  Vernon  H.,  "Design  of  Damping  and  Control  Matrices  by  Modification  of  Eigen¬ 
values  and  Eigenvectors,""  accepted  for  publicatin  in  Journal  of  Shock  and  Vibration, 
submitted  Oct.,  1993. 

PI.  Porter,  B.  and  R.  Crossley,  Modal  Control,  Barnes  and  Noble  Books,  New  York,  1972. 


-92- 


Rl.  Rogers,  L.  C,  "Derivatives  of  Eigenvalues  and  Eigenvectors,"  AIAA  Journal,  8y  943-944, 
(1970). 

51.  Silverberg,  L.,  "Uniform  Damping  Control  of  Spacecraft,"  7.  Guidance.  Control, 
Dynamics,  9(2),  221-227(1986). 

52.  Simoniam,  S.  S.,  "Synthesis  of  Discrete  Passive  Vibration  Dampers,"  Society  of  Auto¬ 
motive  Engineers,  Paper  851935,  Warrendale,  PA.  1986. 

53.  Smith,  David,  "Optimization  of  Damping  in  a  3-Layered  Beam  Using  the  Finite  Element 
Method,"  Senior  Honors  Thesis,  The  Pennsylvania  State  University,  May,  1991. 

54.  Srinathkumar,  S.,  "Eigenvalue/Eigenvector  A.s.signment  Using  Output  Feedback,"  IEEE 
Transactions  on  Automatic  Control,  AC-23,  no  I,  pp.  79-81,  1978. 

55.  Starek,  L.,  and  Inman,  D.  J.,  1991,  "On  Inverse  Problems  with  Rigid  Body  Modes,"  ASME 
Journal  of  Applied  Mechanics,  Vol.  50,  pp.  1 101-1 104. 

56.  Starek,  L.,  Inman,  D.  J.,  and  Kress,  A.,  1992,  "A  Symmetric  Inverse  Vibration  Problem," 
ASME  Journal  of  Vibration  and  Acoustics,  Vol.  1 14,  pp.  564-568. 

VI.  Venkayya,  V.  B.,  and  V.  A.  Tischler,  "Frequency  Control  and  Its  Effect  on  the  Dynamic 
Response  of  Flexible  Structures,"  AIAA  Journal,  23(  1 1 ),  1 768- 1 774  ( 1 985). 

Wl.  Wittrick,  W.  H.  "Rates  of  Change  of  Eigenvalues,  with  Reference  to  Buckling  and 
Vibration  Problems,"  Proceedings  of  the  Royal  Society,  66,  590-591  (Sept.  1962). 


♦ 


-93- 


